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I.  INTRODUCTION 


This  document  is  the  final  report  prepared  under  Contract 
N00014-70-C-0161 ,  Modification  P00005.  The  purpose  of  the  program 
performed  under  this  contract  is  to  continue  the  research  on  10.6-pm 
image  upconversion  detectors  and  to  determine  the  feasibility  of  develop¬ 
ing  a  10. 6-um  range-gated  upconversion  imaging  laser  radar  receiver 
with  150  by  150  resolution  elements.  To  meet  the  program  objectives, 
the  following  areas  were  studied: 

•  Experimental  determination  of  the  effect  of  a 
diverging  pump  on  image  resolution  for  the  image 
space  configuration  upconverter  system  and  com¬ 
parison  of  the  results  with  analytical  prediction 

•  Analytical  determination  of  the  effect  of  bire¬ 
fringence  and  a  truncated  pump  beam  on  image 
resolution 

•  Development  of  techniques  to  handle  power  dissi¬ 
pation  inside  the  nonlinear  crystal 

•  Investigation  of  techniques  for  improvement  of 
optical  efficiency 

An  experimental  model  image  upconverter  capable  of  real  time 
direct  viewing  has  been  assembled  and  tested.  The  optical  arrangement  of 
the  model  was  arranged  to  be  operated  in  either  the  image  space  configura¬ 
tion  (Type  I)  or  in  the  Fourier  space  configuration  (Type  III).  The  experi¬ 
mental  model  can  also  be  operated  either  C\V  or  pulse  gated.  Imaging 
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characteristics  such  as  resolution,  field  of  view,  and  the  upconvr  rsion 
efficiency  of  the  experimental  model  were  measured  for  different  types 
of  operating  modes  and  optical  configurations. 

Analytical  studies  of  the  imaging  characteristics  of  the  para¬ 
metric  image  upconverter  were  also  carried  out  in  the  program.  Effects 
of  crystal  birefringency  and  pump  beam  intensity  distribution  on  the  imag¬ 
ing  properties  were  considered.  A  number  of  useful  relationships  for  sys¬ 
tem  design  was  obtained.  The  experimental  results  and  the  analytical 
results  were  carefully  compared  and  the  excellent  agreement  obtained 
made  it  possible  to  choose  and  design  the  optimum  image  upconverter 
system  for  a  particular  application. 

Some  of  the  material  properties  of  Ag^AsS^  such  as  absorption  co¬ 
efficient  at  1.06  pm  and  the  temperature  distribution  in  the  crystal  under 
laser  beam  irradiation  were  also  studied  experimentally  and  analytically. 

An  objective  of  this  final  report  is  to  compile  and  summarize  infor¬ 
mation  which  is  useful  for  the  understanding  and  design  of  a  parametric 
image  upconversion  system.  A  descriptive  summary  of  the  important 
characteristics  of  the  parametric  image  upconverter  is  presented  in  the 
next  section  for  those  who  are  interested  in  a  synopsis  of  the  results  of 
this  investigation. 


In  Section  III,  system  design  considerations,  design  formulas, 
and  procedures  for  two  types  of  upconverter  systems  are  presented.  Quan 
titative  comparisons  of  the  imaging  characteristics  of  different  types  of 
optical  systems  are  also  contained  in  this  section  and  should  be  sufficient 
for  the  performance  evaluation  of  the  upconverter  system. 

Conclusions  and  recommendations  for  future  work  are  pre¬ 
sented  in  the  last  section. 


II.  SUMMARY  OF  IMPORTANT  CHARACTERISTICS  OF 
PARAMETRIC  IMAGE  UPCONVERTERS 

The  principal  physical  process  responsible  for  parametric 
image  upconversion  is  the  second-order  nonlinear  susceptibility  of  a 
single  crystal  nonlinear  material  such  as  proustite.  The  incoming  object 
waves  and  the  intense  pump  waves  interact  with  each  other  in  the  nonlinear 
material  and  generate  the  sum-frequency  image  waves.  In  order  to  obtain 
efficient  interaction,  the  object  waves  and  the  pump  waves  must  satisfy 
the  phase-match  condition  ii:  the  material.  Among  several  phase-match 
techniques  useful  for  optical  parametric  upconversion,  the  simplest  tech¬ 
nique  suitable  for  imaging  application  is  birefringent  phase-matching  in  a 
single  crystal  bulk  material  since  the  spatial  information  of  incoming 
object  waves  can  easily  be  preserved  with  this  technique.  To  obtain  bire¬ 
fringent  phase-match  in  a  bulk  material,  the  object  waves,  the  pump  waves, 
and  the  upconverted  image  waves  must  be  properly  polarized  and  the  non¬ 
linear  single  crystal  must  be  properly  oriented  with  respect  to  the  system 
optical  axis.  The  phase-match  requirement  and  the  birefringent  nhase- 
match  technique,  a  a  will  be  discussed  later,  give  rise  to  the  imaging 
characteristics  particular  to  the  parametric  image  upconverter  and  re¬ 
stricts  the  parametric  image  upconverter  to  certain  limited  applications. 
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The  phase-match  requirement  and  the  crystal  dispersion  im¬ 
poses  two  basic  constraints  on  the  upconversion  process,  namely: 

•  Usable  spectral  bandwidth  of  the  object  radiation 

•  AUcwable  angular  separation  between  the  object 
waves  and  the  pump  waves 

The  spectral  acceptance  bandwidth  of  parametric  upconverter 
depends  on  the  material  dispersion.  It  is  inversely  proportional  to  the 
interaction  length  and  is  usually  narrow  unless  the  material  dispersion  at 
IR  wavelength  and  sum-frequency  wavelength  is  identical.  Spectral  accept¬ 
ance  bandwidth  of  a  1-cm  long  proustite  pumped  by  a  1.  06- pm  beam  for 
upconverting  at  10.6  um  is  about  6  cm  x.  This  narrow  spectral  acceptance 
bandwidth  of  the  upconverter  using  proustite  virtually  eliminates  using  the 
proustite  image  upconverter  as  a  passive  thermal  image  detector  in  the 
8  to  14  pm  spectrum.  The  application  of  the  proustite  image  upconverter 
system  is  presently  limited  to  the  detection  of  objects  either  radiating 
sharp  spectral  lines  or  illuminated  by  a  laser  beam. 

Spectral  bandwidth  may  be  increased  slightly  by  using  a  con¬ 
verging  or  diverging  pump  beam  since  different  pump  angular  spectrum 
will  phase-match  with  a  different  IR  wavelength.  A  multicolor  thermal 
image  upconversion  might  be  possible  with  this  technique  of  bandwidth 
widening. 
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Because  of  narrow  spectral  acceptance  bandwidth  of  presently 
available  proustite  upconverters,  only  the  image  upconversion  of  mono¬ 
chromatic  objects  is  considered  in  this  report. 

The  acceptance  angle,  defined  as  the  maximum  incidence  angle 
allowed  for  IR  object  waves  with  respect  to  the  phase-match  direction,  is 
inversely  proportional  to  the  square  root  of  interaction  length.  The  pro¬ 
portionality  constant  depends  on  material  refraction  indices  and  three  wave¬ 
lengths  involved.  For  upconversion  of  10.6  um  with  a  1.06-/im  pumped 
1-cm  long  crystal,  the  acceptance  angle  was  about  9  degrees  when  non- 
critically  phase -matched.  The  acceptance  angle  can  be  widened  slightly 
by  orienting  the  nonlinear  crystal  to  phase-match  at  two  different  angles. 

To  restrict  the  angle  of  incidence  of  IR  object  waves  at  the 
crystal  surface  within  the  acceptance  angle  limit,  input  wave  to  the  non¬ 
linear  amterial  must  have  an  f-number  greater  than  the  inverse  of  the 
acceptance  angle  expressed  in  radians.  This  means  an  f-number  greater 
than  6.  5  for  a  1-cm  long  proustite. 

Two  types  of  basic  optical  configurations  suitable  for  para¬ 
metric  image  upconversion  are  shown  in  the  next  section.  In  the  image 
space  configuration  (Type  I)  system,  a  distant  object  (at  IR  wavelength) 
is  imaged  in  the  nonlinear  crystal  by  the  objective  telescope.  The  upcon¬ 
version  process  takes  place  in  the  real  image  space.  An  upcnnverted 


virtual  image  is  relayed  to  the  image  recording  surface  by  the  relay  lens. 
The  upconversion  process  in  the  Fourier  space  configuration  (Type  in) 
takes  place  in  the  angular -frequency  space.  A  real  IR  image  of  the  distant 
object  is  formed  by  the  objective  telescope  in  the  front  focal  plane  of  the 
first  Fourier  transform  lens  in  front  of  the  nonlinear  crystal.  The  object’s 
spatial  Fourier  components  are  upconverted  in  the  nonlinear  material  and 
the  inverse  Fourier  is  transformed  by  the  second  Fourier  transform  lens. 
An  upconverted  real  image  of  the  object  is  then  formed  in  the  back  focal 
plane  of  the  second  Fourier  transform  lens.  In  the  Type  I  system,  IR 
object  rays  originating  from  a  point  interact  with  a  localized  portion  of  the 
pump  beam.  In  the  Type  III  system,  IR  object  rays  from  a  point  interact 
with  the  whole  pump  beam. 

The  acceptance  angle  constraint  on  the  object  waves  introduces 
limitations  on  both  the  field  of  view  (FOV)  and  the  resolving  power  of  the 
parametric  image  upconverter  system  since  the  constraint  affects  the 
spatial  frequency  response  of  the  image  upconverter. 

For  the  Fourier  space  configuration  (Type  III)  system,  the 
acceptance  angle  constraint  restricts  the  size  of  IR  image  field  in  the 
front  focal  plane  of  the  first  Fourier  lens  and  consequently  the  system 
FOV.  The  maximum  size  of  the  IR  image  field  formed  by  the  objective 
telescope  is  equal  to  the  product  of  the  acceptance  angle  and  the  focal 
length  of  the  first  Fourier  transform  lens.  The  system  FOV  is  then  equal 


to  ?,  quotient  of  the  IR  image  field  size  and  the  effective  focal  length  of 
objective  telescope.  The  intensity  distribution  of  upconverted  image  field 
of  a  uniform  intensity  IR  object  field  is  independent  of  the  pump  beam  inten¬ 
sity  distribution  but  is  determined  by  the  phase-match  condition  and  follows 
a  (sin  X/X)  function  where  X  is  a  normalized  transverse  distance  from  the 
optical  axis. 

For  the  image  space  configuration  (Type  I)  system,  the  phase- 
ma*ch  constraint  introduces  thickness  aberration  and  thus  limits  the  image 
resolution.  Linear  resolvable  upconverted  image  line  width  in  the  crystal 
due  to  thickness  aberration  is  proportional  to  the  square  root  of  interaction 
length.  The  proportionality  constant  depends  on  refractive  indices  and  the 
wavelengths.  For  a  1-cm  long  proustite  upconverting  a  10.6-pm  object  by 
the  1.06  uin  pump  beam,  it  is  about  100  um.  Since  the  upconversion  trans¬ 
verse  magnification  of  Type  I  system  is  unity,  the  linear  resolvable  line 
width  cf  an  IR  object  focused  on  the  crystal  is  equal  to  the  upconverted 
image  resolution  width.  This  means  that  the  Type  I  system  would  be 
diffraction-limited  if  the  f-number  of  the  IR  objective  telescope  is  larger 
than  7.7.  When  the  f-number  of  the  IR  objective  telescope  is  matched  to 
the  upconversion  acceptance  angle,  that  is,  the  incidence  angles  of  IR 
object  rays  are  restricted  to  the  acceptance  angle  limit,  the  IR  angular 
resolution  is  independent  of  the  interaction  length  and  is  1.3  times  larger 
than  the  diffraction  resolution  of  the  telescope  aperture. 
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The  effects  on  the  imaging  characteristics  by  nonplanar  pump 
waves,  introduced  either  by  truncation  of  the  plane  wave  pump  by  a  finite 
size  crystal  aperture  or  by  a  nonuniform  transverse  intensity  distribution 
of  the  pump  beam,  are  theoretically  and  experimentally  investigated.  The 
finite  interaction  aperture  due  to  nonplanar  pump  wave  affects  differently 
the  imaging  properties  of  two  types  of  upconverter  system. 

For  the  Fourier  space  configuration  (Type  HI)  system,  the  dif¬ 
fraction  by  the  nonplanar  interaction  aperture  limits  the  system  resolving 
power  if  other  apertures  are  large  compared  to  the  interaction  aperture. 

The  resolving  power  of  a  plane  wave  pumped  infinite  aperture 
in  a  Type  I  system  is  unlimited.  The  IR  Rayleigh  angular  resolution  at 
the  nonlinear  crystal  is  limited  by  the  diffraction  angle  of  the  interaction 
aperture.  If  the  interaction  aperture  is,  for  example,  a  uniform  ampli¬ 
tude  truncated  plane  wave,  the  IR  resolution  angle  is  identical  to  the  dif¬ 
fraction  limit  1.  22  X^/a  where  a  is  the  diameter  of  a  uniform  pump  beam. 
When  the  crystal  diameter  is  large  compared  to  the  pump  beam  size  and 
the  pump  beam  intensity  distribution  is  a  Gaussian,  the  IR  resolution  angle 
becomes  2x^/rra  where  a  is  the  pump  beam  diameter  (1/e  diameter).  It 
must  be  noted  that  the  resolving  power  of  the  Type  III  system  is  very  sensi 
tive  to  the  pump  beam  characteristics.  Use  of  a  pump  laser  radiating  in 
the  lowest  transverse  mode  is,  therefore,  required  for  a  high  resolution 
Type  HI  system. 


The  interaction  aperture  characterized  by  the  intensity  distri¬ 
bution  and  size  of  the  pump  beam  influences  the  intensity  distribution  and 
size  of  the  upconverted  image  field  of  a  Type  I  system.  If  the  object  field 
intensity  distribution  is  uniform,  the  urconverted  image  field  reproduces 
exactly  the  intensity  distribution  of  the  pump  beam.  This  can  easily  be 
visualized  if  one  realizes  that,  in  a  Type  I  system,  the  object  rays  origi¬ 
nated  from  a  point  interacts  with  localized  pump  field.  The  resolving 
power  of  Type  1  system  is  independent  of  the  interacting  aperture  charac¬ 
teristics.  This  implies  that  the  resolution  of  the  Type  I  system  is  insen¬ 
sitive  to  pump  beam  angular  divergence  and  a  multimode  pump  beam  can  be 
used  without  degrading  the  resolving  power  as  long  as  the  angular  diver¬ 
gence  in  the  pump  beam  angular  divergence  is  within  the  pump  acceptance 
angle  limit.  Use  of  a  multitransverse  mode  pump  beam  does  not  degrade 
the  resolution  of  a  Type  I  system.  It  must,  however,  be  avoided  if  pos¬ 
sible,  because  hot  spots  in  the  multimode  pump  beam  might  damage  the 
nonlinear  material  and  voids  in  the  multimode  pump  beam  create  blind 
spots  in  the  image  field.  The  FOV  of  the  Type  I  system  is  independent  of 
the  interaction  length  and  is  equal  to  the  quotient  of  the  effective  inter¬ 
action  aperture  and  the  effective  focal  length  of  the  1R  objective  telescope. 

Conversion  efficiency  (defined  as  a  ratio  of  upconverted  image 
radiant  power  density  and  IR  object  radiant  power  density)  of  parametric 
image  upconverter,  is  not  a  constant  over  the  upconverted  image  field  but 
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is  a  function  of  transverse  position.  For  a  Type  I  system,  the  conversion 
efficiency  is  linearly  proportional  to  the  local  pump  power  density.  For  a 
Type  in  system,  it  is  proportional  to  the  peak  pump  power  density  and  has 
a  (sin  X/X)  functional  dependence  over  the  upconverted  image  field  where 
X  is  a  normalized  transverse  position  of  the  image  field.  The  dependence 
of  upconversion  efficiency  of  the  image  upconverter  on  the  nonlinear  sus¬ 
ceptibility  of  the  material,  the  interaction  length,  the  wavelengths,  and  the 
material  refractive  indices  is  identical  to  that  of  single  resolution  element 
upconverter.  To  maximize  conversion  efficiency  for  a  given  material, 
large  pump  power  density  and  longer  interaction  length  are  required.  Un¬ 
like  the  single  resolution  element  upconverter  that  is  judged  by  the  con¬ 
version  efficiency  alone,  the  figure  of  merit  of  the  multielements  image 
upconverter  should  be  expressed  as  a  product  of  peak  conversion  efficiency 
and  the  number  of  resolution  elements.  The  figure  of  merit  is  linearly 
proportional  to  the  peak  pump  power  and  the  interaction  length  but  is  inde¬ 
pendent  of  the  interaction  aperture  area. 

For  a  1-cm  long  proustite  upconverting  10.6  pm  radiation  by 
interacting  with  a  1.06-pm  pump,  the  normalized  peak  conversion  effi- 
ciency  is  about  5  x  10  per  W/cm  of  pump  beam  density  excluding  other 
optical  losses.  Because  of  low  conversion  efficiency  per  unit  pump-power 
density,  together  with  its  low  damage  threshold,  the  10. 6- pm  proustite 
image  upconverter  system  will  find  its  useful  application  in  range¬ 
gated  pulsed-lidar  imaging  operation.  In  this  particular  mode  of  operation, 
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the  proustite  can  be  Irradiated  bv  a  peak  pump  density  of  1-MvV/cm  without 
causing  material  damage.  The  maximum  peak  upconvcrsion  efficiency  of 
unity  is,  therefore,  theoretically  expected  from  this  proustite  system.  In 
practice,  because  of  optical  loss  and  other  factors  yet  unexplained,  maxi¬ 
mum  peak  efficiency  of  about  10  percent  is  practically  realizable  with 
1  MW/cm  pump-power  density.  Peak  conversion  efficiency  (including 
all  optical  losses)  greater  than  10  percent  is  difficult  with  the  proustite 
because  of  the  material  damage.  Techniques  for  increasing  proustite 
damage  threshold,  such  as  antireflection  coatings,  forced  air  convect  on 
cooling,  and  immulsion  in  index  matched  liquid  coolant  were  found  to  be 
ineffective  under  repetitive  pulsed  pump  laser  irradiation. 

A  slight  deviation  from  the  phase-match  condition  due  to  ma¬ 
terial  temperature  rise  was  observed  under  high  level  CW  laser  irradia¬ 
tion  (30  ~  50  W/cm2).  This  detun  .»g  effect  can  easily  be  corrected  by 
initial  angular  offset.  Angular  detuning  under  pulsed  laser  Irradiation 
was  not  studied  experimentally.  Analytical  consideration  of  heat  transfer 
in  the  material  is  included  in  the  appendixes. 


2-0 


III.  SYSTEM  DESIGN  CONSIDERATIONS 


System  design  of  the  parametric  image  upconvertcr  is  consid¬ 
ered  in  this  section.  Two  basic  optical  configurations  are  considered  and 
a  specific  design  example  of  10.6-um  image  upconversion  is  presented. 
Design  formulas  are  given  without  proof  in  this  section  since  their  deriva¬ 
tion  is  given  in  the  appendixes.  The  basic  conditions  imposed  on  the  up- 
converter  are: 

•  Nonlinear  material  is  placed  in  the  near  field  of  a 
TEM^  mode  pump  beam 

•  Fine  structure  patterns  in  the  upconverted  image 
due  to  the  crystal  birefringence  are  negligible 

•  Noncritical  phase-match  direction  coincides  with 
the  system  optical  axis 

•  Diffraction  due  to  optical  components  other  than 
the  interaction  aperture  is  neglected 

The  first  condition  ib  introduced  since  it  gives  the  best  resolu¬ 
tion  for  the  Type  III  system. 

The  second  condition  is  usuall)  satisfied  for  10.6-pm  upcon- 
verters  using  proustite  because  the  fine  structure  pattern  due  to  birefrin¬ 
gence  is  masked  by  other  dominant  effects  such  as  lens  resolution  and 
image  recorder  resolution. 

The  third  condition  is  a  natural  and  the  most  simple  cnoicr  of 
the  reference  axes. 
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The  fourth  condition  is  imposed  to  isolate  the  image  upron- 
version  characteristics  from  the  properties  of  other  optical  components. 
Inclusion  of  imaging  properties  of  other  optical  components  can  be  easily 
carried  out  by  the  standard  method.  It  is  therefore  rot  considered  in  this 
report. 

A.  UPCONVERTED  IMAGE  FIELD  SIZE  AND  SYSTEM  FIELD  OF  VIEW 

The  size  of  the  upconvcrted  image  field  and  the  system  field  of 
view  ara  governed  by  different  factors  for  the  two  optical  configurations. 

The  phase-match  condition,  consequently  the  IR  acceptance  angle  coo Jnt , 
determines  the  system  FOV  of  the  Typo  Ill  (Fourier  space  configuration)  sys¬ 
tem.  The  interaction  aperture,  determined  by  the  pump  beam  size  or  the 
nonlinear  material  cross  section,  is  responsible  for  lim.iing  the  size  ci 
the  upconvcrted  image  field. 

For  the  Type  I  (image  space  configuration)  system  (Figure  3-1), 
the  size  at  the  IR  object  field  on  the  crystal  is  equal  to  the  crystal  transverse 
dimension  if  the  pump  beam  size  is  considerably  larger  than  the  crystal.  If 
the  cryrtal  is  larger  than  the  pump  beam  size,  the  IR  object  field  size  is 
equal  to  the  pump  beam  size,  taking  the  l/s^  intensity  radius  as  the 
boundary.  The  system  FOV  is  obtained  by  dividing  the  field  size  with  the 


3-2 


IR  OBJECTIVE 
LEWS 

fo 


NONLINEAR 

CRYSTAL 


RELAY  LENS  FOR 
UPCONVERTED  IMAGE 


h 


i_ 

{2 


A.  IMAGE  SPACE  CONFIGURATION 


UPCONVERTED 

IMAGE 


V73-457 


IR  OBJECTIVE 
LENS 


fo 


IR 

I  “AGE 


FIRST 
FOURIER 
TRANSFORM 
LENS  NONLINEAR 
h  CRYSTAL 

on 

fi  — j — fi - 


SECOND 

FOURIER 

TRANSFORM  LENS  Uf>COfr/£RTED 
f2  l*AG£ 


B.  FOURIER  SPAtt  CONFIGURATION 


FIGURE  3-1.  OPTICAL  SYSTEMS  OF  PARAMETRIC  IMAGE  UP- 
CONVERTER 
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effective  focal  length  of  the  objective  telescope.  The  design  formulas 
are 


lil 


If  b  w 


If  b  >  w 


(3-1) 


•o  *  <3-2' 

ahere 

atj  ■  radius  of  !R  object  field  or.  the  crystal 
b  ■  radius  at  the  nonlinear  crystal 
w  ■  pump  beam  radius  (l/e*  Intensity  radius) 

r 

f  ■  effective  focal  length  of  the  objective  telescope 

The  upconverted  Image  field  radius  relayed  by  the  relay  lens  onto  the 
Image  recording  surface  Is  then 

\l  "  1  *11 

uhcre  is  the  transverse  magnification  of  the  relay  lens  (M(  ■  Sj). 
For  Type  111  (Fourier  space  configuration)  system,  the  lit  object  field 
radius  is  given  by 


<1 V2 


(3-3) 
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f j  ■  focal  length  of  1st  Fourier  transfer  lens 
©a  ■  acceptance  angle  of  the  upconverter  process 

The  system  FOV  is  again 


(3-4) 


The  acceptance  angle  is  determined  by  the  phase-match  condition  and  is 
given  by 


2  ViV» 

**np  * 


1/2 


(3-5) 


uhere  n's  are  indices  of  refraction,  x's  are  wavelengths,  and  <  is  the 
interaction  length.  The  radius  cf  upconvcrtcd  image  field  o\cr  the  image 
detecting  surface  is 


2v,v 

MV 


1/2 
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Combining  equations  3-3  and  3-0  th*  system  transverse  mapiification  be¬ 


comes 


B.  RESOLUTION 


The  factor  limiting  the  resolving  power  c i  Type  I  system  is 
thickness  aberration  introduced  by  the  phase-match  constraint.  The  re¬ 
solvable  IR  object  line  width  on  the  nonlinear  crystal  is  independent  of  the 
pump-power  distribution  and  the  interaction  aperture  and  is  given  by 


rl 


VsVj 

Tnsni*p  / 


1/2 


The  system  IR  angular  resolution  Is  then 


(3-8) 


rA 


np  XsM 

- — T 

ffn  n.  X  f 
8  i  p  o 


1/2 


(3-9) 


For  the  Typo  III  system,  the  resolution  ia  limited  by  diffraction  of  the 
interaction  aperture.  If  the  TEM  pnmp  beam  size  is  considerably  larger 
than  the  crystal  diameter  (w  <  b),  the  resolvable  IR  object  line  width  and 

r 

the  angular  resolution  arc  expressed  by 


fj-Xjlj/.Wp  (3-10) 

>’3‘Xl,l/**p,o  ‘3-U> 
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Ii  the  TEM  pump  beam  size  is  considerably  larger  than  the  crystal 


oo 


diameter  tb  «  w  )  and  if  pump  beam  intensity  over  the  crystal  is  fairly 
P 

uniform,  the  IR  object  resolvable  line  width  and  the  angular  resolution 
are  expressed  by 


r3=0. 61  Xjlj/b 
y  3  =  0. 61  X^/b^ 


(3-12) 


(3-13) 


The  number  of  resolution  elements  obtainable  for  an  imaging  system  is 
equal  to  the  square  of  a  quotient  of  the  system  FOV  and  the  angular  resolu¬ 
tion. 


For  the  Type  I  system,  it  is  given  by  combining  equations  3-9 


and  3-2. 


where 


Lil 


N, 


4A11  ns  nl  XP 

1  Vsxi4 


A  ■  irb  (crystal  cross  sectional  area)  if  b  <  <  w 


A?  *  n  Wp4  (pump  beam  area)  is  b  >  wp 


(3-14) 


(3-lo) 
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For  Tyne  HI  system  it  is  obtained  from  equations  3-11,  3-13,  and  3-3 


where 


Nttt  =  - 

III  n 


8A.„  n  n.  X 
i3  s  l  p 


pVi4 


if  b  <  w 

P 


if  b  >  w 
-  P 


(3-16) 


(3-17) 


Comparison  of  equations  3-14  and  3-16  indicates  that  the  resolving  power, 
expressed  in  the  number  of  obtainable  resolution  elements  of  Type  HI  sys¬ 
tem,  is  better  by  a  factor  of  two  than  that  of  the  Type  I  system  if  the  iden¬ 
tical  pump  beam  is  used  for  upconversion. 


C.  CONVERSION  EFFICIENCY 

Conversion  efficiency  of  the  parametric  image  upconverter  re¬ 
quires  a  different  definition  from  that  of  the  single  resolution  element  up¬ 
converter.  For  the  single  resolution  element  upconverter,  it  is  expressed 
by  the  ratio  of  total  upconverted  power  to  total  IR  input  power.  The  image 
upconverter  contains  many  resolution  elements  and  a  nonuniform  intensity 
distribution  within  the  upconverted  image  field.  Conversion  efficiency  of 
the  image  upconverter  must,  therefore,  be  defined  for  a  pair  of  single 
resolution  conjugate  elements.  Since  the  intensity  distribution  (envelope) 
of  the  upconverted  image  field  takes  a  certain  known  functional  form,  only 
the  conversion  efficiency  at  the  center  of  the  upconverted  image  field  is 
required  to  be  considered. 
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There  are  two  possible  ways  of  defining  the  conversion  effi¬ 
ciency.  The  first  definition  is  to  express  it  as  a  ratio  of  upconverted  power 
to  IR  input  power  contained  in  a  pair  of  on-axis  conjugate  resolution  ele¬ 
ments.  This  definiu  n  is  analogous  to  that  of  the  single  resolution  element 
upconverter.  The  second  definition  is  to  express  it  as  a  ratio  of  peak  power 
density  of  a  resolvable  on-axis  upconverted  image  to  peak  power  density  of 
the  conjugate  on-axis  IR  object.  The  latter  definition  is  adopted  in  this 
report  since  it  is  better  suited  to  the  imaging  application. 


The  on-axis  conversion  efficiency  of  two  basic  optical  configu¬ 
rations  is  identical  and  is  given  by: 


W  512  tr5  tdjL  t-2  P 
.  -  s  _  eff  p 

V  ~  nr  IT - - - 2“  X~ 

ir  10  Cnjnsnpxs  p 


(3-18) 


(b  >  w  ) 
-  P 


where 


t 


2 

peak  power  density  in  V  /cm  of  on-axis  upconverted 
image  of  the  conjugate  iR  object  spot  of  minimum 
resolvable  size 

2 

peak  power  density,  in  W/cm  of  on-axis  IR  object 


2  2  2 

(=t.  t  t  )  resultant  optical  surface  transmission  of 

lop 

the  nonlinear  material 


deff  =  effective  SHG  coefficient  along  the  phase-match 
direction  in  cgs  unit 


l  =  interaction  length  in  cm 


pump  beam  power  in  watts 
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=  {l  =  i,  p,  s)  index  of  refraction  along  the  phase-match 
direction  of  nonlinear  material  at  wavelength  denoted 
by  the  subscript 


X  =  sum-frequency  wavelength  in  cm 
s 


2 

=  pump  beam  cross-sectional  area  in  cm 
=  speed  of  light  in  cm/s 


Equation  3-18  is  valid  for  the  case  where  the  crystal  diameter 
is  larger  or  at  least  equal  to  the  pump  beam  diameter.  When  the  pump 
beam  is  considerably  larger  than  the  crystal  diameter,  the  conversion 
efficiency  is  expressed  in  the  same  functional  form  as  that  of  equation  3-14 
but  with  a  different  proportionality  constant.  A  practical  upconverter  sys¬ 
tem  will  satisfy  the  condition  w  <  b,  since  this  condition  is  required  for 

r 

efficient  use  of  the  available  pump  power  and  for  better  resolution.  It  must 
be  noted  that  equation  3-14  is  identical  to  the  conversion  efficiency  expres¬ 
sion  of  the  phase-matched  parametric  upconversion  of  three  plane  waves. 
This  is  no  surprise  if  one  realizes  that  the  condition  along  the  optical  axis 
of  the  image  upconverter  is  almost  the  same  as  that  of  the  phase-matched 
plane -wave  interaction. 


The  on-axis  photon  conversion  efficiency  is,  from  equation  3-14, 
expressed  by 


512  tt5  t  d2  ,  l2  P 
.  a  _  eff  p 

%  "  v  xs/xi  -  77^7 - 


10  '  C  n.  n  n  X  X-  A 

l  p  ssi  p 


(3-19) 
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The  relative  intensity  distribution  of  upconverted  image  field  with  a  uniform 
IR  object  field  emittance  is  considered  next. 


For  the  Type  I  system,  the  upconverted  image  field  envelope 
is  the  same  as  that  of  pump  beam  intensity  distribution  if  the  crystal  is 
larger  or  at  least  equal  to  the  pump  beam  size.  The  size  of  the  upcon¬ 
verted  beam  is  identical  to  the  pump  beam  size.  This  relation  between  the 
pump  beam  and  the  upconverted  image  field  intensity  envelope  holds  for 
non  TEMqo  mode  pump  beam.  Hot  spots  and  voids  in  a  multimode  pump 
beam  will,  therefore,  be  reproduced  in  the  upconverted  image  field  and 
appear  as  blemishes. 


The  upconverter  image  field  intensity  envelope  of  Type  HI  sys¬ 
tem  is,  contrary  to  Type  I  system,  independent  of  the  pump  beam  intensity 
distribution.  It  is  determined  by  the  phase-match  condition  and  is  ex¬ 
pressed  by 


WsW 
ws  (p8  =  o)  - 


2  n  n.  X  X  f„ 
s  ips  2 


\ 


I 


(3-20) 


where 

pg  =  radial  distance  in  the  upconverted  image  field 
from  the  optical  axis 

W  (p i  )  =  power  density  of  upconverted  image  field  at  p 

S  S  c 

f„  =  focal  length  of  second  Fourier  transform  lens 
(Figure  3- IB) 
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Equation  3-6,  which  describes  the  size  of  the  upconverted  image  field,  is 

obtained  from  equation  3-20  by  taking  the  first  root  of  the  sine  function  as 

the  limiting  size  ana  solving  for  the  radial  distance  o  . 

s 

The  conversion  efficiency  achievable  with  proustite  is  about  a 
factor  of  less  than  the  theoretical  calculation  expressed  by  equation  3-14. 
Reasons  for  the  discrepancies  between  experimental  measurements  and 
analytical  predictions  are  not  known  with  certainty.  Possible  reasons  are 
crystal  absorption  and  scattering.  Some  discussion  on  absorption  and 
scattering  in  proustite  is  presented  in  Appendix  R. 

D.  SYSTEM  SIGNAL- TO-NOISE  RATIO 

To  achieve  any  advantage  from  upconversion,  detection  or 
recording  of  the  upconverted  image  with  a  reasonable  quantum  efficiency 
in  a  signal-shot-noise  limited  mode  is  desirable.  Presently  available 
image  detectors  suitable  for  the  image  upconverter  application  are  either 
S-20  or  S-l  photocathode  image  intensifiers.  An  S-20  photocathode  is 
better  than  an  S-l  response  cathode,  but  its  use  requires  a  Ruby  laser 
pump.  To  achieve  a  signal-shot-noise  limited  operation,  the  image  in- 
tensifier  tube  must  be  either  gated  or  cooled  or  both  to  reduce  dark  noise. 
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The  narrow  spectral  acceptance  of  the  upconverter  is  an  ad¬ 
vantage  when  the  upconverter  system  is  used  as  a  laser  receiver,  since  it 
lessens  the  background  noise.  Parametric  uptonversion  does  not  (to  a 
first  order  approximation)  generate  internal  noise.  Noise  at  the  sum- 
frequency  originates  from  the  background  thermal  radiation  within  the  IR 
acceptance  spectral  bandwidth. 

Other  sources  of  noise,  which  must  be  added  to  the  background 
noise  are  the  pump  leakage  noise  due  to  imperfect  pump  suppression  and 
the  image  detector  dark  current  noise. 

The  signal-to-noise  ratio,  which  depends  on  the  radial  distance 
from  the  optic  axis,  of  the  image  upconverter  system  using  an  image  intcn- 
sifier  tube  for  direct  viewing  is  then  expressed  by: 


(3-21) 

where 

rjd  =  quantum  efficiency  of  the  image  intensifier  photocathode 

rj(p)  ■  conversion  efficiency  given  by  equation  3-14  multiplied 
by  a  radial  distance  dependent  function 

Wj  *  on-axis  radiant  omittance  of  IR  object 
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t  "transmittance  at  Xr  of  output  optical  components 
t jQ  ■  transmittance  at  x  4  of  input  optical  components 

Wb  *  thermal  background  radiant  omittance  at  the  input 
surface  of  the  nonlinear  material 

Oj  ■  collection  efficiency  of  1R  telescope 

W  (o)  ■  pump  beam  power  density  (depends  on  radial  distance) 

r 

t  *  transmittance  at  of  output  optical  components 

\Vt  ■  noise  equivalent  irrad  Lance  of  the  photocathode  of 
image  intensifier  tube 

t  ■  observation  time  (integration  time) 

f  ■  pulse  repetition  rate 

Tj  ■  pulse  w  idth  of  active  IR  illumin.  ior 

r  ■  width  of  receiver  gate 
B 

It  Is  assumed  in  equation  3-21  that  the  width  of  the  receiver  gate 
Is  equal  to  the  pump-pulse  width  and  the  pump-pulse  width  Is  longer  than  the 
IR  illuminator  pulse  width.  If  the  pump  pulse  is  shorter  than  the  IR  pulse 
width,  t i  in  equation  3-21  should  be  replaced  by  t  .  For  a  system  oper- 

*  r» 

ating  in  C\V  mode,  the  conditions  r  ■  t  and  fr  "1  hold. 

*  i  g  g 

To  illustrate  the  relative  magnitude  of  noises,  a  proust ite  up- 

converter  pumped  by  a  TEMq0  mode  Nd:YAG  laser  beam  Is  considered 

below.  It  was  assumed  that  the  Image  detector  Is  a  cooled  S-l  photo- 

cathode  image  intensifier  tube.  Noise  equivalent  input  energy  W,  of  typical 
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cooled  S-l  tube  is  approximately  10" 11  W/cm2  at  0.967  wm.  The  thermal 
background  emittance  of  the  upconverter  system  is  generated  within  the 
nonlinear  crystal  if  it  is  not  cooled.  If  the  input  IR  optical  components  and 
the  crystal  were  cooled,  the  thermal  background  would  originate  in  the 
background  scene. 


where 


2  h  &Vi  cc 

(e  1  C 


-1) 


v h  j 

atio  hy.  At. 

(e  1  b-l) 


for  room  tempera¬ 
ture  nonlinear  crystal  (a) 

(3-22) 

for  cooled  IR  optics 
and  crystal  (b) 


Tc  -  crystal  temperature 
Tfa  =  background  temperature 


r  =  crystal  emissivity  (»0.05  for  proustite  at  10  um) 
c 

For  a  room  temperature  proustite  with  bv^/C  -  0.05  cm-1,  wb  « 

•6  2 

4  x  10  W/cm  .  For  t  =  0.25,  the  quantity  representing  background 

oU 

noise  w.t  is  about  10"®  w/cm2. 
b  so  ' 


The  quantity  representing  pump  leakage  noise  in  equation  3-21 
becomes  independent  of  absolute  pump  power  density,  since 


M  /  -2CW  / 

topWp(o)/”(p)  =  2Wpo‘op  e  P/K  <*»>  ‘  “op6  y  «•> 
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where 


g(p)  =  e 


-2p2/w2 


P  for  Type  I 


g(o)  =  sine2 


I  ff<np  xi^‘ 


l  2nsniXpXsf2  I 


for  Type  HI 


Furthermore  g(o)  for  Type  in  system  can  be  expressed  approximately  by 

2/2 
-p  A 

e  r  by  proper  choice  of  since  a  sine  function  can  be  approximated 
by  Gaussian  function  up  to  the  first  zero  of  a  sine  function.  If  this  is  done, 
pump  leakage  noise  contribution  will  be  independent  of  the  radial  distance 
and  becomes 


Wp<‘>)topA’(p)  =  2top/K 


(3-23) 


where  K  is  expressed,  from  equation  3-14 


K  = 


512  „5tc  de„2  l2 


x-7 


10  Cn.n  n  xe 
X  s  p  s 


2~  («5  x  10  for  a  1-cm  proustite) 


The  quantity  representing  photocathode  dark  noise  is 


Wt/n(p)=2Wt/W  Ke 


_  2/  2 
•2p  /wp 


(3-24) 
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The  conditions  that  the  background  noise  will  dominate  the  pump  leakage 
noise  and  the  tube  dark  noise  are  then 

t  „  <  t  W  k/2  (=  2.  5  x  10-13  for  1-cm  proustite  and  t  =  0.25) 
op  so  b  r  so  ' 

and  (3-25) 

o  2/  2 

-2p  /w  ,  2 

We  2  W. /K  t  W.  (=  40  W/cm  for  proustite  and 

p°  t  so  b  cooled  S-l  cathode 

with  t  =0.  25) 
s 

The  above  results  indicate  that,  for  the  system  considered,  a 

2 

minimum  of  126-dB  pump  rejection  and  a  minimum  of  40  W/cm  '  pump- 

2 

power  density  at  1/e  radius  are  required  for  background  or  signal-shot- 
noise  limited  operation.  When  tQp  =  125  dB  and  Wp  =  40  W/cm3,  the  back¬ 
ground  noise,  pump  leakage  noise,  and  the  tube  dark  noise  are  approxi¬ 
mately  equal. 

The  above  condition  for  the  pump  suppression  can  be  satisfied 
by  use  of  two  interference  filters  and  a  polarization  'ilter.  The  trans¬ 
mission  at  x  through  these  filters  would,  however,  not  exceed  0.3  in 
s 

practice.  The  choice  of  tgo  =  0.25  is,  therefore,  appropriate  for  the 
system  design. 
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Insufficient  pump  suppression  in  front  of  the  image  intensifier 
affects  the  system  performance  differently  for  two  different  types  of  optical 
configuration.  For  the  Type  I  system,  the  pump  beam  is  expanded  by  the 
image  relay  lens  and  appears  larger  than  the  original  size  in  the  image 
plane.  Since  the  power  density  of  an  expanded  pump  beam  is  less  than  the 
original  pump  beam,  required  pump  suppression  for  this  system  is  less 
than  above  calculation.  The  leaked  pump  wave  increases  background  noise 
over  entire  image  field  for  Type  I  system.  The  pump  beam  is  focused  on 
a  point  in  the  image  plane  by  the  second  Fourier  transform  lens  for  the. 
Type  HI  system.  Insufficient  pump  suppression  would,  therefore,  create 
a  blind  spot  in  the  image  field  where  the  leaked  pump  beam  is  focused. 
Required  pump  suppression  for  this  system  is  more  stringent  than  Type  I 
system  if  one  wishes  to  eliminate  a  blind  spot  within  the  FOV. 


Since  the  background  noise  dominant  conditions  can  be  satis¬ 
fied  for  a  pulse-gated  upconverter  system,  the  image  upconverter  is  either 
background  noise  limited  or  signal- shot-noise  limited.  For  this  case  the 
cn-axis  signal-to-noise  ratio  becomes 
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(3-26) 


where  R_ 
po 

2 

joule/cm  . 
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is  the  on-axis  pump  energy  density  per  pulse  in 
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It  is  assumed  in  equation  3-26  that  the  illuminator  pulse  width 


is  equal  or  longer  than  the  pump-pulse  width. 


E.  A  SYSTEM  DESIGN  EXAMPLE  FOR  10.6-ym  IMAGE  UPCON- 
VERTER  USING  PROUSTITE _ 

A  design  of  a  pulse-gated  active  image  upccnverter  using  a 
Nd:YAG  laser  pumped  proustite  is  considered.  The  following  assumptions 
are  made: 


•  A  cooled  S-l  image  intensifier  tube  is  used  for  image 
viewing 

•  Pump  leakage  and  image  tube  dark  noise  are  small 
compared  to  the  thermal  background  noise 

•  Proustite  can  be  safely  irradiated  with  peak  power 

o 

density  of  0.  5  MW/cm  with  repetitively  pulsed 
pump  beam 

•  Image  degradation  by  the  optical  components  is 
negligible 

•  All  lenses  are  AR  coated  so  that  reflection  losses 
from  lenses  are  negligibly  small 

•  Phase-match  change  due  to  crystal  temperature 
rise  is  corrected  if  there  is  any 

•  Pump  beam  is  in  TEM  mode 

oo 

•  Gate  open  duration  is  equal  to  or  shorter  than  the 
illuminator  pulse  width  (t.  =  t  ) 

•  D 

•  The  crystal  diameter  is  larger  or  at  least  equal  to 
the  pump  beam  diameter 

•  Proustite  is  at  room  temperature 
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Parameter  values  used  In  the  following  design  calculation  arc 
listed  as  follows: 
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The  required  pump  beam  sire  and  peak  pump  power  arc  ealeu 
la  ted  first  from  equations  3-14,  3-20,  and  assumed  peak  pump  on-axis 
power  density  of  0. 5  MW /cm2. 
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Using  parameter  values  given  previously,  one  obtains 


-3  */2 

wpI  ■=  5.S3  x  10  J  (Nt  ) 

-3  */2 

Wpm'3.91  x  10  d(Nf) 

where  subscripts  I  and  III  denote  Type  I  and  Type  111  systems,  respectively. 
For  N  *  150  x  150  and  i  *  1  cm 

w  j  =  0.829  cm 
and 


WpUj  c  0.  587  cm 

The  peak  pump  power  required  to  give  0.  5  MW/cm2  cn-axis  peak  power 
density  is  then  0. 54  MW  for  the  Type  I  system  and  0. 27  MW  for  the  Type  HI 
system. 

The  theoretical  on-axis  peak  conversion  efficiency  is,  from 
equation  3-14,  equal  to  0. 125  for  0.  5  MW/cm2  on-axis  pump  power  den¬ 
sity.  Practically  achievable  conversion  efficiency  is  about  one-half  of  the 

theoretical  estimate,  therefore,  practical  on-axis  conversion  efficiency 
2 

with  0.  5  MW/cm  pump  density  is  about  6  percent. 


Figure  3-2  is  a  plot  01  equation  3-26  for  a  pulse-gated  10.6-um 
proustite  system  with  parameter  values  given.  For  direct  viewing,  the 
integration  time  t  is  approximately  equal  to  0.2  s,  the  integration  time 
of  human  eyes.  To  avoid  flicker,  the  pulse  repetition  rate  should  be  at 
i.east  16  PPS.  Pump  pulse-width  of  1/2  microsecond  to  1  microsecond 
should  be  used  for  range-gated  operation  since  the  range  gate  of  75  to 
150  m  is  desirable  for  practical  radar  application.  The  quantity  RpfT 
for  the  system  considered  ranges  from  0.8  to  1.6.  The  on-axis  minimum 
discernible  object  power  density  is  about  10“ 5  W/cm2  for  unity  contrast 
objects. 

The  off-axis  minimum  discernible  object  power  density  is  ob¬ 
tained  by  multiplying  the  on-axis  value  by  a  factor  equal  to  inverse  of 
amplitude  distribution  function  of  the  upconverted  image  field.  For  a 
Type  I  system,  the  image  field  amplitude  follows  a  Gaussian  distribution 
and  therefore  the  minimum  discernible  object  power  density  at  the  edge 
of  the  image  field  is  e  times  of  on-axis  value,  that  is,  6. 4  x  10"b  W/cm2. 
For  a  Type  III  system,  the  amplitude  distribution  of  image  field  is  ex¬ 
pressed  by  a  (sin  X/X)  function.  Minimum  discernible  object  power  den¬ 
sity  at  the  transverse  position  corresponding  to  the  root  of  sin  X/X  func¬ 
tion  becomes,  of  course,  infinity. 
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FIGURE  3-2.  PLOT  FOR  A  PULSE-GATED  10.6-uM  PROUSTITE  SYSTEM 


! 


The  power  density  scattered  back  from  a  laser  illuminated 


target  is 


P  r 
t  t 


W.  =  1  1 

l  o?  R2 


(3-27) 


where 


Pt  =  laser  illuminator  peak  power 
rt  =  target  reflectivity 
0^  =  illuminating  beam  divergence 
R  =  range 

Ideally  the  illuminator  beam  divergence  is  matched  to  the  re¬ 
ceiver  field  of  view.  The  required  illuminator  power  per  unit  range  can 

be  calculated  from  equations  3-26  and  3-27  for  a  given  signal-to-noise 

-5  2 

ratio.  For  0^.  =  0. 1  rad,  r^.  =  0. 1,  andW^  =  2.72  x  10  .  W/cm  ,  the 

4  2 

minimum  illuminator  peak  power  per  1  km  range  is  2.1  x  10  W/km  . 

The  illuminator  energy  per  pulse  required  for  a  5-km  range  is  then  0.  5  J, 
which  can  be  generated  by  a  reasonable  C02  laser. 

Performance  characteristics  and  optical  parameters  of  the 
receiver  system  considered  are  summarized  in  Table  3-1. 
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TABLE  3-1.  CHARACTERISTICS  OF  PROUSTITE  10.6-uM  RANGE¬ 
GATED  RADAR  IMAGE  UPCONVERTER  SYSTEM 


r 

Optical  System 

I 

III 

U 

Nonlinear  material 

1-cm  long  proustite 

1-cm  long  proustite 

0 

FOV 

5  degrees  to  1/e^ 

5  degrees  null  to  null 

Number  of  resolution  elements 

150  x  150 

150  x  150 

u 

Frame  rate 

16  ft/s 

16  ft/s 

r 

Image  detector 

Cooled  S-l  inten- 
sifier 

Cooled  S-l  inten- 
sifier 

L 

Object  wavelength 

10. 6  jim 

10.6  um 

i" 

Upc  on  verted  wavelength 

0.967  um 

0.967  um 

1; 

Pump  wavelength 

1 . 06  u m 

1 . 06  u  m 

L 

On-axis  pump  power  density 

0.  5  MW/cm2 

0.  5  MW/cm^ 

On-axis  conversion  efficiency 

0.125 

0.125 

h 

Required  minimum  pump 
beam  diameter 

1. 7  cm 

1.  2  cm 

Peak  pump  power 

0.54  MW 

0.27  MW 

r» 

Gate  width 

1  JIS 

1  u  s 

Pump  energy  per  pulse 

0.  54  J 

0.27  J 

r. 

On-axis  minimum  discernible 
object  power  density 

10" 5  W/cm2 

10" 5  W/cm2 

i. 

Objective  lens  effective 
F-number 

7 

7 

i: 

System  optical  efficiency  at 
sum  frequency 

0.25 

0.25 

Range 

5  km 

5  km 

il 

COg  illuminator  laser  energy 

0.  5  J 

0.  5  J 

per  pulse 

-* 

Target  reflectance 

0.1 

0.1 

it  » 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

Experimental  and  analytical  studies  carried  out  in  this  program 
have  demonstrated  that  the  proustite  image  upconverter  pumped  by  a 
1.06  fim  pump  beam  is  capable  of  reproducing  an  image  of  10.6  um 
illuminated  objects  with  resolution  exceeding  150  x  150  elements.  It  has 
also  demonstrated  that  a  range-gated  radar  image  upconversicn  with  a 
reasonable  range  can  be  constructed  with  presently  available  components. 

The  imaging  characteristics  and  parameters  restricting  resolu¬ 
tion  and  efficiency  are  identified.  A  set  of  design  formulas  useful  for  the 
system  design  is  obtained.  The  appropriate  applications  of  image  upcon¬ 
verter  techniques  at  present  are: 

•  Range-gated  radar  image  system 

•  IR  holography 

•  Laser  diagnostic  technique 

•  Laser  weapon  aim-point  monitoring 

In  order  to  advance  the  state  of  the  art  of  image  upconversion 
technique  development,  the  following  technical  areas  are  recommended: 

•  Development  of  new  nonlinear  material  with  improved 

characteristics 
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Development  erf  high  quantum  efficiency  low  noise 
image  tube  with  a  narrow  spectral  response  at  x 

Development  of  a  high  efficiency  band-pass  filter 
with  sufficiently  large  out-of-band  rejection 

Improvement  of  pump  laser  characteristics  and 
laser  construction. 
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APPENDIX  A 


REAL  TIME  IMAGE  UPCONVERSION  EXPERIMENTS  WITH  THE 
DIRECT  VIEWING  LABORATORY  MODEL 

An  experimental  model  upconverter  capable  of  real  time  direct 
viewing  of  upconverted  images  has  been  set  up  during  this  period. 

The  imaging  characteristics  of  the  upconverter  were  measured 
first  with  the  proustite  crystal  pumped  by  a  CW  TEMqo  mode  beam.  The 
resolving  power  of  this  model  arranged  in  the  Fourier  space  configuration 
was  measured  to  assess  how  close  it  came  to  achieving  diffraction  limited 
operation.  The  model  was  then  set  up  in  the  image  space  configuration  to 
measure  the  extent  of  thickness  aberration. 

The  measured  angular  resolution  of  approximately  1  milliradian 
for  the  Fourier  space  configuration  was  within  a  factor  of  two  of  the  Rayleigh 
resolution  of  a  uniformly  illuminated  object  lens  aperture.  The  number  of 
resolution  eltments  contained  in  a  10-degree  field  ot  view,  which  is  also 
experimentally  checked,  is  consequently  better  than  the  resolution  objective 
for  the  program  of  150  by  150  elements.  Upconverted  images  of  a  10.  6-ym 
illuminated  USAF  1951  resolution  chart  containing  groups  smaller  than  1  lp/.nm 
were  obtained.  Because  the  resolution  chart  does  not  fill  the  field  of  view, 
the  actual  number  of  resolvable  elements  contained  in  the  recorded  images 
which  fills  a  7.  5-degree  angle  is  114  by  114. 
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The  resolving  power  of  the  model  arranged  in  the  image  space 
configuration  was  also  measured.  The  experimental  results  prove  (as 
expected)  that  the  resolvable  image  linewidth  is  limited  by  the  thickness 
aberration.  The  measured  image  resolution  linewidth  of  60  4  m  agrees 
with  the  physical  optics  prediction  of  68  pm,  but  is  approximately  three  to 
four  times  narrower  than  the  geometric  optics  calculation.  The  main  reason 
for  failure  of  the  geometric  optics  theory  to  accurately  predict  the  extent  of 
aberration  is  probably  due  to  the  inability  of  the  geometric  optics  to  account 
for  the  intensity  variation  of  the  upconverted  image  rays.  Upconverted  images 
of  the  resolution  chart  with  62  by  62  resolution  elements  were  obtained  with 
the  experimental  model  in  the  image  space  configuration.  Measured  angular 
resolution  is  slightly  worse  than  the  diffraction  limited  angular  resolution  of 
the  object  lens. 

The  effect  of  diverging  pump  beam  on  the  resolving  power  of  the 
upconverter  system  in  the  image  space  configuration  was  investigated.  The 
fact  that  resolution  of  the  image  space  configuration  system  using  a  multi  - 
mode  pump  beam  does  not  degrade  from  that  of  a  TEMoq  mode  pumpec'.  sys¬ 
tem  was  verified  experimentally.  However,  a  nonuniform  intensity  corre¬ 
sponding  to  the  multimode  pattern  is  superimposed  on  the  upconverted  image. 

linage  upconversion  experiments  with  the  image  space  configured 
model,  using  pulsed  pump  beam  and  a  gated  image  tube,  were  also  carried 
out.  Some  preliminary  results  of  the  experiments  are  also  presented. 
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1. 


EXPERIMENTAL  ARRANGEMENT 


The  experimental  model  upconverter  system  is  shown  schematically 
in  Figure  A-l.  The  system  was  arranged  so  that  the  optical  configuration  and 
the  pump  beam  parameters  could  be  readily  changed  with  minimum  physical 
alteration  of  the  setup.  A  2. 5-cm  diameter,  and  0. 6  cm  thick  proustite  crys¬ 
tal  was  used.  The  upconverter  was  aligned  ior  Type  II  noncritical  phase 
matching.  The  proustite  was  oriented  such  that  its  crystallographic  X  axis 
lies  in  the  plane  of  paper  and  is  parallel  to  the  front  surface.  Its  optic  axis 
lies  in  a  plane  perpendicular  to  the  paper  and  makes  an  angle  of  20  degrees 
with  the  normal  (along  the  C  axis)  to  the  polished  front  and  back  surfaces. 

A  continuously  pumped  Nd:YAG  laser  radiating  electric  fields,  polarized 
parallel  to  the  plane  of  paper,  was  used  as  a  pump  souree.  The  laser  is 
fitted  with  an  acousto -optic  Q-switch  allowing  pulsed  as  well  as  CW  opera¬ 
tion.  The  beam  parameters  of  the  1.06-pm  pump  beam  emanating  out  of  the 
beam  expander  are  shown  in  Table  A-l  for  various  modes  of  operation. 

TABLE  A-l.  1.06  Mm  PUMP  BEAM  PARAMETERS 
AFTER  BEAM  EXPANSION 


MODE 

AVERAGE 

POWER 

BEAM 

CROSS  SECTION 

BEAM 

DIAMETER 

BEAM 

DIVERGENCE 

CW  TEMco 

5.0  watts 

Gaussian 

1. 1  cm 

0. 16  mrad 

CW  Multi- 
mode 

7.  5  watts 

1. 1  cm 

1 . 0  mrad 

Q-Switched 

TEMoo 

1000  PPS 

0.  5  watt 

Gaussian 

1. 1  cm 

0. 16  mrad 
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FIGURE  A-l.  SCHEMATIC  DIAGRAM,  EXPERIMENTAL  MODE 
IMAGE  UPCONVERTER  SYSTEM 


A  polarizer  external  to  the  Nd:  YAG  laser  removes  field  com¬ 
ponents  polarized  along  the  y  direction.  The  chopper  reduces  the  average 
power  of  the  pump  beam  by  serving  as  an  attenuator  when  required.  Glass 
lenses  Lj  and  collimate  and  expand  the  pump  beam. 

A  multimode  CC>2  laser  radiating  with  an  electric  field  polarized 
perpendicular  to  the  plane  of  the  paper  was  used  to  back  illuminate  the  trans¬ 
parent  1951  USAF  resolution  chart.  The  C02  beam  diameter  at  the  resolu¬ 
tion  chart  is  sufficiently  large  to  illuminate  the  entire  resolution  chart  by 
adjusting  the  10. 6-pm  path  length  by  means  of  mirrors  as  shown  in  Figure  A1 
The  position  of  BaF2  lens  Lg,  relative  to  the  proustite  crystal  and  the  object, 
is  determined  by  the  configuration  of  the  system  (image  space  or  Fourier 
space). 

The  pump  wave  and  the  object  wave  are  combined  by  a  dichroic 
mirror.  An  improved  dichroic  mirror  was  obtained  during  this  reporting 
period.  This  mirror  reflects  more  than  99  percent  of  the  YAG  beam  and  has 
a  transmission  loss  at  10. 6  pm  of  about  13  percent.  As  shown  in  Figure  A-l, 
the  angle  of  incidence  of  the  YAG  beam  on  the  dichroic  mirror  is  about  20°. 

This  eliminates  potential  vignetting  of  10.6-pm  object  waves  by  the  mirror 
holder . 

A  Gian  air  polarization  filter  and  narrowband  interference  filters 
are  used  for  suppression  of  1.06-pm  pump  wave.  Transmission  through 
these  filters  at  0.967  pm  is  about  10  percent.  Lens  L 
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is  used  to  form 


upconverted  real  images  on  the  photo  surface  of  the  image  intensifier  tube. 
The  image  intensifier  tube  is  an  ITT  Model  4740  three-stage  gated  tube  with 
an  S-l  photocathode.  Its  pertinent  characteristics  are  listed  in  Table  A-2. 
The  image  intensifier  tube  can  be  jooled  to  -30°C  with  dry  ice  to  reduce 
tube  thermal  noise. 

TABLE  A-2.  CHARACTERISTICS  OF  ITT  MODEL  F4740  IMAGE  TUBE 


Useful  Cathode  Diameter 
Radiant  Sensitivity  at  0. 967  pm 
Center  Resolution 
Peripheral  Resolution 
Paraxial  Magnification 
Distortion 
Luminous  Gain 

Conversion  Efficiency  at  1.06  pm 

Radiant  Power  Gain  Peak  Wavelength 

Radiant  Power  Gain  1 . 06  p 

Equivalent  Background  Energy  Input  at 
Peak  Wavelength 

Equivalent  Background  Energv  Input  at 

1.06p 


23  mm 
1.3  inA/W 
28  lp/mm 

25  lp/mm  at  7. 0-ram  radius 
0.84 

18  percent  at  10-mm  radius 
6000  fL/fc 

'■°00-000  wTTSgf 

1040  W/W 
208  W/W 
2  x  10-9  W  7cm2 

8  x  IQ'9  W'cm2 


To  check  the  resolving  power  of  the  0.967  pm  image  detector 
used  in  the  system,  consisting  of  lens  L^,  the  interference  filter,  the  image 
intensifier  tube,  and  the  camera,  the  resolution  chart  is  illuminated  with 
visible  light  and  imaged  on  the  photosurface  of  the  image  intensifier  i«'be. 
Figure  A-2  shows  a  typical  image  recorded  on  the  film.  The  smallest 
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RESOLUTION  LIMITS  OF  THE  0.9fc?-;im  IMAGE 
DETECTOR  OF  THE  EXPERIMENTAL  MODEL  IS  37 
AS  SEEN  IN  THE  ILLUSTRATION 


FIGURE  A-2.  IMAGE  OF  THE  RESOI.l  ITON  CHART  ILLUMINATED 
BY  VISIBLE  LIGHT 
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resolvable  linewidth  on  the  output  phosphor  of  the  image  tube  is  37  pm. 

A  modified  Polaroid  scope  camera  was  used  to  record  the  upconverted  images 
produced  by  the  image  tube. 

2.  IMAGE  UPCONVERSION  IN  FOURIER  SPACE 

The  limiting  resolving  power  of  the  experimental  model  was 
measured  by  operating  the  model  in  Fourier  Space  configuration  and  pump¬ 
ing  the  nonlinear  material  with  a  TEMq()  mode  beam.  Diffraction  limited 
image  upconversion  is  possible  for  this  configuration  if  the  pump  is  a  plane 
wave.  Thickness  aberrations  are  negligible  for  this  configuration. 

A  1951  USAF  resolution  chart  was  placed  in  the  front  focal  plane 
of  lens  L3.  The  distance  between  BaF2  lens  Lg  and  the  proustite  crystal  is 
120  mm  which  is  equal  to  the  focal  length  of  lens  Lg.  Lens  L  is  placed  so 
that  the  proustite  crystal  and  the  photocathode  of  the  image  intensifier  tube 
coincide  with  the  front  and  back  focal  planes  of  lens  L 4,  respectively.  The 
focal  length  of  lens  L4  is  580  mm. 

The  pump  beam  diameter  at  the  proustite  front  face  is  1. 1  cm  and  its 
beam  divergence  is  0. 16  rnrad.  The  pump  beam  is  in  TEMqo  transverse  mode. 
Pictures  reproduced  from  some  of  the  upconverted  images  and  recorded 

by  the  Polaroid  camera  are  shown  in  Figure  A-3.  The  smallest  element  that 
can  be  visually  resolved  in  these  figures  is  the  group  2/l.  Figure  A-3  A  has 
better  resolution  along  the  x  direction,  but  Figure  A-3B  has  better  resolution 
along  the  y  direction.  (See  Figure  A-l  for  the  coordinate  system. )  Inspection 
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FIGURE  A-3. 


of  upconvcrted  images  on  the  output  phosphor  surface  of  the  image  intensifier 
tube  indicates  that  group  1/6  is  resolved  along  both  the  x  and  the  y  directions. 
The  system  angular  resolution  is  calculated  by  the  formula 

y  =  hir/f3  (A-l) 

where 

hir  =  resolvable  linewidth  of  the  test  chart 
f3  =  effective  focal  length  of  lens  L„ 

O 

The  measured  angular  resolution  corresponding  to  the  1/6  and  2/1  resolu¬ 
tion  groups  is : 

Group  1/6  h.r  =  1/7. 13  mm  =1.17  x  10'3  rad 
Group  2/1  hir  =  1/8  mm  =  1.04  x  10"3  rad 

Measured  angular  resolution  of  1. 17  mrad  is  2.  5  times  larger  than  the 
Rayleigh  resolution  angle  of  3-cm  diameter  lens  L,. 

Resolution  degradation  from  the  Rayleigh  resolution  angle  of  the 
object  lens  is  mainly  due  to  nonuniform  transverse  intensity  distribution  of 
the  pump  beam.  It  will  be  shown  later  that  a  Gaussian  profile  of  the  pump 
beam  would  introduce  an  additional  Gaussian  attenuation  factor  on  the  angular 
spectrum  of  the  upconverted  image  field.  It  is,  therefore,  expected  that  the 
resolving  power  of  the  experimental  model  is  limited  by  the  Gaussian  window 
diffraction  rather  than  the  object  lens  diffraction. 
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The  system  field  of  view  was  measured  by  translating  the  object 
along  the  x  and  y  directions  until  the  edge  of  the  test  chart  disappeared. 

The  measured  system  FOV  of  0. 183  rad  (10.  5°)  agrees  with  the  previous 
result.  Figure  A-3A  shows  that  the  FOV  covers  more  than  the  area 
occupied  by  the.  resolution  chart.  The  number  of  resolvable  elements  con¬ 
tained  in  the  FOV  is  found  to  be  156  by  156,  if  one  takes  group  1/6  as  the 
minimum  resolvable  element.  The  image  shown  in  Figure  A-3A  contains 
approximately  114  by  114  resolution  elements. 

The  measured  transverse  magnification  of  0.99  agrees  within 
10  percent  of  the  estimated  value  of  1. 08,  which  is  calculated  by  the  ex¬ 
pression 
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where 
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uc 


As  f4 
S  -^=0.44 


Xir  f3 


:  upconversion  magnification 


(=2.94):  camera  magnification 
M^,  (  =  0.  81):  image  tube  magnification 

3.  IMAGE  UPCONVERSION  IN  IMAGE  SPACE 


Image  upconversion  experiments  were  carried  out  with  the 
experimental  model  arranged  in  the  image  space  configuration.  Lens  L„  of 

w 

Figure  A-l  was  positioned  such  that  a  real  10.6-um  image  of  object  is  focused 
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onto  the  crystal.  The  photosurface  of  the  image  intensifier  and  the  proustite 
crystal  coincide  with  the  conjugate  planes  of  lens  L4> 

It  will  be  shown  analytically  in  Appendix  C  that  the  resolvable  image 
linewidth  of  the  image  space  configuration  upconverter  pumped  with  a  plane 
wave  is  limited  by  thickness  aberration.  It  is  the  purpose  of  this  experiment  to 


measure  the  aberrations  and  to  determine  the  limiting  resolution  of  the  image 


space  configured  upconverter.  The  Nd: YAG  laser  was  operated  in  the  TEM 

oo 


mode  to  minimize  the  effect  of  diverging  pump  beam. 


Figure  4  shows  reproduced  pictures  of  upconverted  images  re¬ 
corded  by  the  camera.  The  magnification  of  image  forming  lenses  and 
are,  respectively,  0.304  and  2.  76.  The  smallest  resolvable  group  (1/3)  has 
a  linewidth  of  2.  52  lp/mm.  The  corresponding  image  width,  calculated  by 
dividing  the  resolvable  target  width  by  the  magnification  factor  of  lens  L,,  is 

O 

60. 3  jim. 


An  experimentally  measured  resolvable  image  linewidth  of  60  jim 
is  approximately  a  factor  of  3  narrower  than  the  resolvable  image  width  pre¬ 
dicted  by  the  geometric  optics  calculation.  As  will  be  explained  later,  the 
physical  optics  calculation  predicts  that  the  minimum  resolvable  image 
width  for  the  given  parameter  values  is  68  jim.  The  experimental  result 
agrees  better  with  the  physical  optics  prediction  than  with  the  geometric 
optics  calculation. 
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In  the  geometric  optics  model,  the  analysis  is  based  on  the 
assumption  that  each  infinitesimal  slab  of  material  is  capable  of  forming  an 
independent  upconverted  image.  The  concept  of  thickness  aberrations  is, 
therefore,  introduced  by  considering  the  upconverted  image  as  a  congregation 
of  images  formed  by  infinitesimal  slabs  of  upconverter  material.  It  has  been 
previously  shown  (by  the  Fourier  optics  treatment)  that  the  aberration  is 
physically  introduced  by  imperfect  phase  matching  of  angular  spectra  of  the 
interacting  electric  fields.  The  angular  spectrum  with  higher  spatial  fre¬ 
quency  suffers  higher  attenuation  due  to  the  phase  mismatch.  The  degree  of 
phase  mismatch  is  proportional  to  the  material  thickness.  In  the  physical 
optics  model  the  upconverted  image  is  formed  by  linear  superposition  of  the 
upconverted  rays  (angular  spectra)  emanating  from  the  exit  surface  of  the 
material.  It  is  expected  that  the  aberrations  calculated  by  these  two  models 
should  be  different  from  each  other.  The  experimental  results  should  agree 
better  with  the  Fourier  optics  predictions  rather  than  the  geometric  optics 
because  the  model  is  in  better  agreement  with  the  actual  physical  process. 

The  number  of  resolution  elements  contained  in  the  image  shown 
in  Figure  A-4  is  approximately  62  by  62  which  is  approximately  a  factor  of  two 
less  than  that  obtained  by  the  Fourier  space  configuration  system.  Measured 

system  magnification  is  approximately  10  percent  greater  than  the  theoretical 
estimate. 

The  system  angular  resolution  can  be  calculated  by  dividing  the 
measured  resolvable  image  width  by  the  distance  between  the  crystal  and 
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FIGURE  A-4. 


UPCON VERTED  IMAGE  OF  1951  USAF  TEST  CHART 
PRODUCED  HY  IMAGE  SPACE  UPCONVERTER  WITH 
SINGLE  MODE  PUMP 
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lens  Lg.  Measured  system  angular  resolution  of  0.4  mrad  is  very  close  to 
the  Rayleigh  resolution  angle  of  0. 42  mrad  of  lens  Lg.  The  experimental 
model  is,  therefore,  nearly  diffraction  limited. 

It  must  be  noted  that  the  Gaussian  distribution  of  the  pump  beam 
does  not  degrade  the  resolution  as  expected  from  the  physical  optics  analysis. 

It  can  easily  be  shown  that  the  Gaussian  amplitude  factor  of  the  angular 
spectrum  of  the  upconverted  image  fields  becomes  unity  for  all  angular 
frequencies  in  the  image  space  configuration  where  the  apparent  IR  object 
location  coincides  with  the  nonlinear  material. 

The  primary  effect  of  Gaussian  profile  of  the  pump  beam  appears 
as  intensity  reduction  of  upconverted  images  toward  the  edge  of  FOV.  The 
size  of  the  10. 6-pm  image  of  the  object  formed  at  the  crystal  would,  there¬ 
fore,  be  considerably  smaller  than  the  pump  beam  size  if  one  wants  to  avoid 
intensity  reduction  at  the  outer  edge  of  FOV.  In  the  experiment,  the  size  of 
the  10. 6-|im  image  of  the  resolution  chart  at  the  crystal  was  3.  6  mm  by  3. 6  mm 
(considerably  smaller  than  the  pump  beam  diameter).  Loss  of  brightness  and 
resolution  at  the  edge  can  even  be  seen  in  Figure  A -4,  which  was  taken  with 
some  precaution  to  minimize  the  Gaussian  window  effect.  If  the  pump  beam 
is  expanded  to  1-inch  diameter  and  usable  pump  power  is  increased  by  four 
times  that  of  the  present  system,  the  number  of  resolution  elements  could 
be  increased  to  a  minimum  of  130  by  130,  since  the  resolvable  image  width 
remains  the  same. 
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It  was  proved  experimentally  that  the  resolvable  image  width 
was  unchanged  while  the  magnification  of  lens  Lg  did  change. 

4.  IMAGE  SPACE  UPCONVERSION  WITH  A  MULTIMODE  PUMP 
BEAM _ _ _ _ _ 

The  purpose  of  this  experiment  is  to  prove  that  a  multimode 
pump  beam  can  be  used  in  the  image  space  configuration  system  without 
sacrificing  the  resolving  power  from  the  single  mode  pumped  image  space 
system.  It  is  widely  recognized  that  the  image  resolution  degrades  with  an 
increase  of  pump  beam  divergence.  This  has  been  demonstrated  experi¬ 
mentally  for  the  Fourier  space  configuration  upconverter  and  for  the 
optical  system  in  which  an  apparent  IR  object  is  located  at  finite  distance 
from  the  nonlinear  crystal. 

The  reduction  of  image  resolution  due  to  the  pump  beam  diver¬ 
gence  (or  equivalently  the  pump  beam  profile)  has  been  analyzed  by  the  ray 
tracing  technique  and  the  physical  optics  approach.  It  has  been  shown  that 
the  amplitude  of  the  angular  spectrum  oi  the  upconverted  image  of  an  on 
axis  point  object  contains  an  attenuation  factor  that  depends  on  the  pump 
beam  intensity  distribution  and  the  object  position.  This  attenuation 
factor  becomes  independent  of  the  pump  beam  profile  when  the  apparent  IR 
object  plane  coincides  with  the  equivalent  pump  beam  plane.  The  resolving 
power  then  becomes  insensitive  to  the  pump  beam  profile  and  will  be  limited 
by  the  thickness  aberration.  It  must  be  understood,  however,  that  the 
intensity  distribution  of  the  upconverted  image  depends  strongly  on  the  pump 
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beam  intensity  distribution.  In  order  to  prevent  undesirable  image  intensity 
variation  across  the  FOV  in  the  image  space  configuration  system,  a  ran¬ 
domly  phased  multimode  pump  beam  is  preferred  over  a  beam  in  a  single 
(or  a  few)  higher  order  modes  which  may  have  holes  in  the  beam. 

Figure  A-5  shows  reproduced  pictures  of  an  upconverted  image  of 
the  resolution  chart.  The  image  shown  in  Figure  A-5A  was  obtained  with  a 
TEMoq  mode  pumped  image  space  system.  Figure  A-5B  shows  an  image  ob¬ 
tained  with  the  same  system  pumped  with  a  multimode  beam.  The  two  images 
remain  essentially  the  same,  indicating  no  degradation  of  resolution  due  to 
the  pump  beam  divergence.  It  was  observed,  however,  that  the  image  ob¬ 
tained  from  the  multimode  pump  beam  was  less  bright  than  that  produced  by 

the  TEM  pumped  beam.  Quantitative  comparison  of  image  brightness  for 
oo 

these  two  cases  was  not  carried  out. 

It  appears  that  dependence  of  upconversion  efficiency  upon  the 
pump  beam  divergence  requires  further  analytical  and  experimental  investiga¬ 
tion  to  design  an  efficient  image  upconversion  system. 

5.  IMAGE  UPCONVERSION  WITH  A  PULSED  PUMP  LASER  AND  A 
SYNCHRONOUSLY  GATED  IMAGE  TUBE _ 

An  experimental  model  pulse-gated  image  upconversion  system 
was  assembled,  and  some  preliminary  experiments  were  carried  out.  The 
optical  layout  of  the  model  is  identical  to  the  image  space  configuration  (Fig¬ 
ure  A-l).  A  block  diagram  of  electronics  circuits  used  for  the  synchronous 
pump  pulse  generation  and  gating  of  the  image  tube  is  shown  in  Figure  A-6. 
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A.  TEM  MODE  PUMP 
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FIGURE  A-5.  UPCONVERTED  IMAGE  OF  1951  USAF  TEST  CHART 
PRODUCED  BY  IMAGE  SPACE  UPCONVERTER 
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Nd:  YAG  LASER 


TO  IMAGE  INTENSIFIEP 
GATE  GRID 


2-4811 


FIGURE  A-6.  BLOCK  DIAGRAM  OF  SYNCHRONOUS  GATE 
ELECTRONICS 


The  Nd:YAG  las  ;r  operating  in  single  transverse  mode  was  Q-switched  by 
an  acousto -optical  Q  switch.  The  pulse  repetition  rate  is  controlled  by  the 
trigger  pulse  generator  and  can  be  varied  from  1  to  5  kHz.  Pump  pulses 
with  200-ns  duration  and  approximately  5-kW  peak  powe»-  were  expanded  to 
1-cm  diameter.  No  damage  to  the  crystal  was  observed.  The  high  voltage 
gate  pulse  generator  has  a  40-ns  rise  time  and  an  80-ns  fall  time.  Pulse  width 
and  the  pulse  position  of  the  gate  pulses  can  be  adjusted  for  the  optimum 
image  signal-to -noise  ratio. 

The  images  obtained  with  the  synchronously  gated  image  tube 
are  shown  in  Figure  A-7A.  The  image  shown  in  Figure  A-7B  is  obtained  by 
the  same  system,  but  the  image  tube  was  not  gated.  The  improvement  in 
the  image  signal-to-noise  ratio  attained  by  gating  the  image  is  quite  remark¬ 
able.  Optinr.zation  of  the  experimental  system  for  range -gated  operation 
was  not  attempted  to  obtain  Figure  A-7A. 
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A.  SYNCHRONOUSLY  GATED  IMAGE  TUBE 


B.  UNGATED  IMAGE  TUBE 


FIGURE  A-7.  UPCON VERTED  IMAGE  WITH  A  PIMP  BEAM 
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APPENDIX  B 

UPCONVE RSION  EFFICIENCY  MEASUREMENTS 
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Conversion  efficiency  of  the  upconversion  process  has  been 

measured  and  compared  to  a  theoretical  calculation  based  on  the  plane  wave 

12  3 

theory  by  several  authors  ’  ’  .  Measured  efficiency,  the  conversion 
efficiency  formula  used  for  checking  their  results,  and  parameter  values 
such  as  SHG  d  coefficients  are  inconsistent  among  investigators.  In  this 
section,  an  attempt  will  be  made  to  clarify  some  of  the  confusion  by  con¬ 
sidering: 

•  Proustite  SHG  d  coefficients 

•  Experimental  measurement  of  upconversion  efficiency 
under  controlled  experimental  conditions 


1.  PROUSTITE  SHG  d  COEFFICIENTS 

The  only  reported  measurement  of  SHG  d^  coefficients  of 
proustite  was  performed  by  K.  H.  Hulme  et  al  by  comparing  SHG  of  1. 152  (im 

4 

with  proustite  to  SHG  with  LiNbO^  along  the  crystallographic  X-axis  .  The 

ratio  Idggl/Msil  of  proustite  was  determined  by  Hulme  et  al  by  the  fre- 

4 

quency  summing  technique  .  Their  measurement  did  not  determine  the 
absolute  signs  of  the  coefficients.  Their  measured  results  are 
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1.6  ±0.1 


(B-l) 
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The  inconsistency  in  proustite  dgi  is  resulted  from  the  spread 
in  LiNbOg  dgl  coefficient  reported  in  the  literature  as  shown  below. 

G.  Boyd5 


d31(LN)  =  (10.6U.0)d36(KDP) 
£ 

R.  Miller  and  Savage0 

d31(LN)  =  (n.9n.7)d36(KDP> 

Bjorkholm  ® 


dJLN>  =  (10.9  ±1.7)  d  (KDP) 


‘31 


36 
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Byer  and  Harris 
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d  (KDP) 
a14 


10.6 


d  (KDP) 
a36 


Careful  review  of  the  above  reported  measurements  indicates  that  Bjorkholm's 
measurement  seems  more  accurate  than  others. 

Using  Jerphagnon  and  Kurtz*9  and  Francois **  measurement 
on  dgg  of  KDP  crystal, 


d  (KDP) 
36 


=  1. 67  x  10~9  esu 
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and  Bjorkholm's  d 


one  finds 


31 


(LN) 


'31 


(P) 


(P) 


*22 


=  28.4  d36(KD)  =  42.6  x  10"9  esu 


=  45.  5  d36(KDP)  =  68.  2  x  10"9  esu 


(B-2) 


Hulme  et  al  gave  estimated  d  coefficients  from  their  measurement  as 


id 


(P) 


22 


=  50  d36(KDP)  =  75  x  10"9 


esu 
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31 


(P) 


=  30  d36(KDP)  =  45  x  10"9  esu 


The  reason  for  this  modification  from  experimental  results  was  not  given  by 
Hulme. 

12 

Boggett  et  al  reported  the  effective  d  coefficient  of  Type  I 
phase  matched  proustite  by  measuring  doubled  10. 6  jum  output  power.  Their 
results  for  0  =22.5°. 


60  d36 


(KDP) 


is  higher  than  the  calculated  value  based  on  Hulme's  estimated  d  coefficients 


d22<P>  *  d36(KDP) 


and  d 


31 


(P)  _ 


=  30  d 


(KDP) 


36 


.  If  one  solves  the  equation 


deff  =  d22  cos  8m  +  d31  sin  8m 
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which  applies  to  Type  I  phase  matching,  assuming  dQ0/d(}1  =1.6  for 

ol 

proustite,  one  finds  d31<P)  =  32.2  d36<KDP)  from  the  Boggett  measurement. 
Boggett's  result  is  very  close  to  the  Hulme’s  estimated  dQ1  value. 

The  values  ^22  =  ^  *  10  csu  and  d^  =  42. 6  x  10  ^  esu  will 

be  used  for  following  system  performance  calculation,  although  Hulme's 
estimated  d  coefficients  are  widely  adapted  by  other  workers. 

The  choice  is  dictated  by  the  conversion  efficiency  measurements 
performed  at  AIL.  The  result  of  the  most  recent  conversion  efficiency  mea¬ 
surement  is  presented  in  paragraph  2. 

2.  CONVERSION  EFFICIENCY  MEASUREMENTS 

The  main  purpose  of  the  experiment  is  to  assess  the  achievable 
conversion  efficiency  of  a  practical  image  upconverter  as  compared  to  the 
theoretical  estimate.  The  experimental  conditions  for  efficiency  measure¬ 
ments  must  be  carefully  controlled  in  order  that  the  experimental  results 
can  be  meaningfully  compared  to  the  theory. 

CW  conversion  efficiency  of  the  experimental  model  was  measured 
with  TEMqo  mode  signal  beam  and  TEMqo  mode  pump  beam.  The  upcon- 
verted  signal  was  detected  by  a  low  noise  SI  photomultiplier  with  accurately 
calibrated  quantum  efficiency.  Optical  efficiency  of  optical  components  such 
as  interference  filter,  optical  windows,  and  the  beam  combining  dichroic 
mirror  were  carefully  calibrated.  The  pump  beam  was  chopped  with  a  duty 
factor  of  one  quarter  to  prevent  excessive  crystal  surface  oxidation. 
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The  input/output  linearity  was  checked  so  that  an  efficiency  mea¬ 
surement  was  taken  within  the  dynamic  runge  of  the  system.  Figure  B-l  shows 
the  plot  of  input  COg  power  versus  the  upconverted  signal  induced  photoelec¬ 
tron  at  the  photocathode.  The  observed  power  upconversion  efficiency  is 
calculated  by 

N  hC 

^u  =  ?T d  P.  X  (B~4) 

p  lr  s 

where 

Ng  =  photoelectron  flux 
r?p  =  quantum  efficiency  of  the  photomultiplier 
h  =  Planck's  Constant 
C  =  velocity  of  light 
Xg  =  wavelength  of  upconverted  signal 
Pir  =  IR  input  power 
d  =  duty  factor  of  light  chopper 
Using  the  following  values: 

d  =  0.25 

“X~=  2-05  *  10  19  Joules  for  \  ,  =  0.967  urn 
s  s 


RELATIVE  C02  POWER 


FIGURE  B-l.  UP (X)N VE RSION  LINEARITY  PHOTON  COUNTS 
VERSUS  RELATIVE  C0o  POWER 


the  measured  system  peak  power  conversion  efficiency  is  4.05  x  10“9  which 

is  approximately  one-third  of  the  theoretical  estimated  efficiency  of 
-8 

1. 46  x  10"  calculated  by  the  expression,  which  applies  to  upconversion  of 
a  TEMqo  1R  beam  with  a  TEM^o  mode  pump  beam, 


Pa  104.44  d  2  l}  P  t  *  <.2P 

„  =  _»  - Si L - E.JJP  .  2. 162  *  10-6°P  ,PP 

'prSl* 


n.  n  n  (A  +  A.) 
i  l  s  p  s  p  i 


(B-5) 


where 

Pg  =  total  upconverted  power  in  watts 

Pj  =  total  IR  power  in  watts 
Pp  =  pump  power  in  watts 

l  =  interaction  length  in  cm 
n^  =  indices  of  refraction  at  (l  =  i,  p,  s) 

Ap  =  pump  beam  area  in  cm^ 

A.  =  IR  beam  area  in  cm2 

t  =  optical  transmission  of  the  experimental  model 

Parameter  values  used  in  the  efficiency  estimate  are: 

deff  =  d22  C0S  °m  =  68  <cos  20°>  *  10"9  esu  =  64  *  10' "9  esu 

n.  =  2.68 

rip  =  2.817 

n  =2.804 
s 

*■=0.6  cm 

Xg  =  0.967  4m 
2 

Ap  =  1. 13  cm  (beam  diameter  =  1.2  cm) 
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2 

=  0. 783  cm  (beam  diameter  =  1  cm) 

P  =  4.8  W 
P 

t  =  7.44  x  10"3 
op 

The  most  probable  cause  for  the  discrepancy  between  the  mea¬ 
sured  efficiency  and  the  prediction  based  on  equation  B-5  is  the  phase  mis¬ 
match  due  to  the  beam  divergence  and  internal  pump  beam  scattering. 
Equation  B-5  was  derived  for  the  case  of  perfectly  phase-matched  inter¬ 
action  of  two  TEMqo  mode  beams.  Other  factors  affecting  efficiency  re¬ 
duction  are  the  beam  walk-off,  the  crystal  birefringence,  and  the  crystal 
absorption. 

CW  system  conversion  efficiency  under  the  experimental  con¬ 
dition  that  nearly  realizes  the  plane  wave  interaction  model  was  measured. 
The  measured  conversion  efficiency  is  approximately  a  factor  of  1.  75  less 
than  the  theoretical  prediction  calculated  by  equation  B-5.  In  this  measure¬ 
ment,  the  10. 6-jim  beam  was  spatially  filtered  by  a  small  pinhole  placed 
at  beam  center.  The  filtered  10.6-ym  wave  front  is  very  close  to  plane  wave 
front.  Beam  divergence  effect  is  minimized  by  aligning  the  optics  so  that  the 
interaction  of  IR  waves  with  the  pump  waves  is  limited  to  the  center  of  the 
pump  beam. 

Conversion  efficiency  of  pulse -gated  upconverter  was  also  mea¬ 
sured  to  ascertain  that  there  is  no  abnormality  in  this  mode  of  operation. 

The  Nd:YAG  laser  was  acousto-optically  Q-switched  with  a  pulse  repetition 
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rate  of  1  kHz.  Pulse  width  and  pulse  shape  of  pump  laser  outputs  were 

measured  with  a  high-speed  Ge  photodetector  of  known  responsivity.  Half- 

power  pulse  width  and  the  peak  power  is  250  ns  and  2.  5  kW,  respectively. 

Upconverted  peak  power  and  the  pulse  shape  are  measured  witn  another  Ge 

photodiode  of  known  responsivity.  Pump  beam  power  density  at  the  center 

2 

of  beam  was  about  5  kW  cm  .  Transmittance  of  optical  components  was 

improved  to  0. 28  for  the  gated  model.  Measured  upconversion  efficiency  of 

-5  -4 

4. 1  x  10  is  about  one-third  of  the  estimated  value  1.1  x  10  calculated 

by  equation  B-5.  The  accuracy  of  above  conversion  efficiency  measurement 

is  about  50  percent  since  the  calibration  of  optical  attenuator  used  for  the 

pump  power  measurement  could  not  be  calibrated  better  than  50-percent 

accuracy.  The  result  obtained  for  pulse  pumped  upconverter  agrees  with 

that  of  CW  upconverter. 


Conversion  efficiency  of  image  upconversion  must  be  defined 
differently  from  equation  B-5  since  there  are  many  resolution  elements  in 
the  upconverted  image.  It  can  be  shown  that  the  conversion  efficiency 
defined  as  a  ratio  of  upconverted  power  density  to  IR  object  power  density 
for  a  pair  of  conjugate  points  is  given  by: 


W  (p  ) 
s  -  s' 


(B-6) 
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where: 


t  -  total  surface  transmittance  at  the  crystal  surfaces 
F(pg)  -  Envelope  factor  which  depends  on  the  optical  configuration 
and  pump  beam  envelope 

Ps  *  radial  distance  of  image  point  measured  from  the  optic 
axis 

*  radial  distance  of  the  conjugate  object  point 

For  a  TEM^  mode  beam  pumped  Type  I  system,  the  envelope 
factor  is  expressed 


(B-7) 


where 

a  =  pump  beam  radius 


For  a  TEMqo  mode  beam  pumped  Type  III  system,  we  have 


n  I  v  t  n  X .  p 

Fm(Ps>  =  Sine2  (  —  P  1  s 


2  n  n  X  xtf  , 
si  p  s  2  / 


where  f ^  =  focal  length  of  second  Fourier  transform  lens. 


(B-8) 
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APPENDIX  C 

THEORETICAL  ANALYSIS 


1.  INTRODUCTION 

Imaging  theory  of  parametric  image  upconversion  was  treated 
previously  by  the  geometric  optics  approach.  An  equation  relating  locations 
of  object,  pump,  and  image  was  derived  and  the  existence  of  image  aberra¬ 
tions  due  to  the  finiteness  of  the  nonlinear  material  thickness  was  made 
clear.  In  geometric  optics  treatment  of  parametric  image  upconversion, 
some  higher  order  factors  that  would  affect  imaging  properties  were  ne¬ 
glected.  It  was  considered  also  in  the  geometric  optics  theory  that  upcon- 
verted  image  is  formed  by  constructive  interference  between  the  waves 
radiated  from  each  differential  slab  of  polarization  in  the  material.  This 
concept  leads  to  thickness  aberration  of  the  image  upconversion.  To  verify 
the  concept  of  constructive  interference  and  to  understand  the  higher  order 
effects,  such  as,  anisotropy,  phase  mismatch,  and  finite  interaction  aper¬ 
ture,  the  problem  must  be  treated  in  the  formalism  of  physical  optics.  It 
is  the  purpose  of  this  study  to  treat  the  problem  in  the  framework  of  scalar 
wave  theory.  The  problem  considered  here  is  to  determine  the  nature  of 
the  image  field  arising  from  a  dielectric  polarization  generated  at  the  sum 
frequency  by  interaction  of  the  object  field  and  the  pump  field  in  the  nonlinear 
material.  The  imaging  properties  of  the  parametric  image  upconversion 
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have  been  treated  by  Flrcster  ^  in  the  formalism  of  Fourier  optics 
that  resembles  closely  the  treatment  of  linearly  invariant  filter  system. 

The  advantage  of  such  an  approach  is  that  the  process  of  upconvcrsion  of 
the  complex  nonplanar  objects  field  to  the  complex  nonplanar  image  field 
can  be  studied  by  Investigating  the  interaction  of  planar  angular  spectrum 
components  in  the  nonlinear  material.  The  results  obtained  by  Fourier 
optics  treatment  describe  more  quantitatively  the  imaging  properties  of 
upconversion  process  and  can  be  tied  directly  to  the  formalism  of  optical 
transfer  function  description  of  an  imaging  system. 

General  formulation  of  image  upconvcrsion  process  taking  into 
account  of  phase  mismatch  and  the  crystal  anisotropy  is  presented  in  para¬ 
graph  2.  The  treatment  presented  follows  Kleinman*1^  closely  and  is  more 

(13) 

rigorous  than  Firester'  The  basic  equations  that  describe  the  angular 
spectrum  of  image  field  in  terms  of  object  and  pump  angular  spectrum  are 
obtained.  The  concept  of  constructive  interference  of  radiation  from  dif¬ 
ferential  slabs  is  also  verified  in  paragraph  2.  Imaging  properties,  such  as 
angular  spectrum  of  image  field,  image  position,  and  the  point  spread  func¬ 
tion  of  two  image  upconversion  systems,  namely  Fourier  space  configuration 
system  and  image  space  configuration  system,  are  considered  in  paragraphs 
3  and  4.  Imaging  with  a  high  order  mode  pump  and  the  effects  of  birefringence 
on  the  image  properties  are  discussed  in  paragraphs  5  and  6. 


C-2 


2.  GENERAL  FORMULATION 
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l. 


a.  BASIC  EQUATIONS 

Mathematically,  the  formulation  is  based  upon  Maxwell's  equa¬ 
tion  for  a  lossless  anisotropic  medium  free  of  charges  and  currents  but  con¬ 
taining  a  prescribed  source  in  the  dielectric  polarization.  All  linearly  induced 
polarization  is  assumed  to  be  included  in  the  dielectric  constant  tensor.  The 
analysis  follows  very  closely  to  Kleinman’s  treatment*15’  l6)  of  second- 
harmonic  generation  except  that  the  analysis  will  be  carried  out  in  Fourier 
spatial  frequency  space  rather  than  in  spatial  coordinate  space.  It  can  be 
shown  that  the  Fourier  decomposed  angular  spectrum  of  the  field  disturbance 
can  be  considered  as  plane  waves  propagating  in  different  directions.  Para¬ 
metric  interaction  of  complex  (nonplanar)  waves  in  the  nonlinear  material 
can  be  understood  by  studying  the  interaction  of  their  angular  spectrum  com¬ 
ponents.  The  Fourier  optics  formalism  puts  the  problem  in  the  framework 
of  linear  transferse-spatial  invariant  system. 


The  theory  of  parametric  upconversion  involving  three  plane 
waves  has  been  treated  in  considerable  detail*17,  18’  19).  under  matching 
conditions  the  plane  wave  signal  field  E.(z)  and  the  output  plane  wave  sum 
frequency  field  E  (z)  satisfy  the  relations* 18 \ 
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E  (z)  = 
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u) 


CO 


1/2 


E.(0)  sin 


it  z 
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n. 
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1/2 


(C-l) 


where  the  characteristic  length  l  for  power  conversion 


2  co  .  co 
1  s 


C  <K,K  ) 


P2  Ep(0) 


;c-3) 


with  K  the  wave  number,  co  the  angular  temporal  frequency,  X  the  appro¬ 
priate  component  of  the  second-order  polarization  tensor,  and  Ep(0)  the 
pump  electric  field  at  input  surface  of  the  nonlinear  material.  In  the  limit 
z  <<  -t,  the  sum-frequency  field  is  given  by^*^, 

1  w 

Es(z)  «  2ux  —  -jr  E.(0)  E  (0)  z  (C-4) 

s  ^ 

(15) 

As  shown  by  Kleinman'  ,  the  linear  growth  of  E  (z)  is  ascribed  to  con- 

s 

structive  interference  between  a  forced  wave  generated  by  the  polarization 

and  a  free  wave  generated  at  the  surface  by  the  forced  wave.  It  is  equally 

valid,  as  we  shall  see  later,  to  ascribe  the  growth  to  constructive  interference 

between  the  waves  radiated  from  each  differential  slab  of  polarization  in  the 

material.  We  developed  geometric  optics  theory  according  to  latter  view. 

Equation  C-4  is  the  solution  for  the  matching  case  satisfying  E  (0)  =  0. 

s 

Actually  equation  C-4  does  not  give  the  correct  sum-frequency  field.  The 


value  of  E  (0)  is  related  to  the  amplitude  of  nongrowing  waves  in  whe 
s 

medium  and  a  reflected  wave  at  the  surface^ These  vaves  are, 
however,  many  orders  of  magnitude  weaker  than  the  growing  wave. 

The  rigorous  analysis  of  upconversion  process  is  based  upon 
the  inhomogeneous  vector  wave  equation^  ®’ 


a' 

L 

CO  | 
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E  (r)  =  4n 
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s 

C  1 

P(r) 

(C-5) 

where  e  =J.+  4tt  ^ ^  is  the  linear  susceptibility  tensor]  and  Eg  are 

respectively  dielectric  constant  tensor  and  electric  field  at  the  frequency 

oo  g  =  to.  +  GOp.  The  plane  polarization  wave  P  (r,  t)  is  related  to  the  signal 

electric  field  £.(r)  and  pump  field  E  (x)  by  the  relation 
*  P 


-ico  t  -i(oo.  +  co  )t 

P(r,  t)  =  P(r)  e  S  =  E  (r)  E  (r)  e  1  P  (C-6) 

—  r 


For  nonplanar  signal  and/or  pump  fields,  which  is  the  case  for 
imaging,  it  would  be  advantageous  to  solve  the  inhomogeneous  wave  equation 
for  the  angular  spectrum  component  waves  rather  than  solving  equation  C-5 
directly.  The  wave  equation  for  Fourier  components  wave  is 


V  xvxfi  (K  ,z)-€  :S  (K  z)  —  ^  if 

—  o  '  —  e7  '  ao  —  c  '-c;  ' 


O) 


£(*,*)  (C-7) 
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It  is  apparent  that  operation  v  on  angular  spectrum  only  has  g  opera 
tion.  In  this  analysis  Cartesian  coordinates  x,  y,  z  are  defined  such  that 
z  is  in  the  propagation  direction  normal  to  the  surface  of  the  crystal,  and 
the  optical  axis  (consider  only  uniaxial  crystals)  is  in  the  xz  plane  (Fig¬ 
ure  C-l).  If  both  signal  and  pump  fields  are  ordinary  wave  polarized  with 
the  electric  vectors  perpendicular  to  the  xz  plane  (along  the  y  direction) 
the  effective  component  of  the  polarization  wave  and  the  sum-frequency 
electric  field  are  polarized  along  the  x  axis.  For  this  simple  case  equa¬ 
tion  C-7  becomes 


1  a2  ' 

la  1 

2 

1  bz2  +  'ej 

Ssx(Ss>z)*in 

s 

c  1 

*xiSS’z) 

(C-8) 

where 

fe  = 

e  e  u 

fx  =  x2iy(Is,2)®6py(«s,2) 

The  general  solution  of  equation  C-7  consists  of  an  inhomogeneous 
wave  (or  forced  wave)  which  is  some  convenient  particular  solution  of  in¬ 
homogeneous  equation,  and  suitable  free  waves  which  are  solutions  of  the 

Homogeneous  equation  obtained  by  setting  jP  =0.  The  free  waves  are 

s 

required  so  as  to  satisfy  the  boundary  conditions  at  the  crystal  surface. 


CRYSTAL  SURFACE 
PARALLEL  TO  x  y 
PLANE 


FIGURE  C-l.  RELATION  BETWEEN  THE  COORDINATES  (x,  y,  z)  AND 
THE  CRYSTALLOGRAPHIC  AXES  (X,  Y,  Z) 


The  free  waves  have  the  form 
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where  the  refractive  index  n(  is  a  function  of  the  free  wave  propagation 
direction  s  (unit  vector)  of  the  wave  vector 

O! 

K  =  -^ngS  (C-10) 

The  relation  between  the  wave  vector  K  of  free  wave  and  the  wave  vector 
Kg  of  the  polarization  wave  P  ( r)  will  be  considered  later.  The  polarization 
direction  Q  satisfies  the  relation 


where 


a  :  *26  =  0 
~s  ~ 
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=s 


( 


n  I  -  0  -  K  k 
s~  -  s  s 


(C-12) 


where  I  is  the  unit  dyadic  and  c.  is  the  dielectric  constant  tensor  at  sum 

frequency  w  .  Equations  C-ll  and  C-12  determine  both  n  and  the  polar- 

s  s 

ization  direction^.  We  define  a  quantity  nr,  called  effective  index,  and 
a  unit  vector  o  by  the  relation 


K 
— s 


CO 

S  /  „ 
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(C-13) 


In  general  n  t  n  (£),  but  n  -  n  (cr)  along  the  phase-match  direction, 
o  s  — 
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When  n'  ^  ng(fl),  the  polarization  wave  is  said  to  be  mis¬ 
matched,  the  mismatch  is  described  quantitatively  by  0  in  the  relation15^ 

n/a-ns  =  0N  (C-14) 

O 

or 

to 

Ks-K=-^0N  (C-15) 

where  N  is  the  direction  (unit  vector)  normal  to  the  incident  surface  of  the 
crystal.  These  equations  determine  both  0  and  s  to  be  associated  with 
a  prescribed  Kg.  Here  ng  and  s  specify  a  free  wave  expressed  by  equation 
C-9  which  must  exist  in  the  presence  of  the  forced  wave  in  order  to  satisfy 
boundary  condition  at  the  crystal  surface.  The  dominant  effect  of  presence 
of  both  forced  and  free  waves  is  the  production  of  growing  wave  through  in¬ 
terference  of  a  free  wave  and  forced  wave^15\  When  0  =  0  the  waves  are 
said  phase  matched  and  the  growing  wave  is  given  by  equation  C-4.  More 
generally  under  nearly  matching  conditions  following  Kleinman^15^  it  can 
be  shewn  that  the  growing  wave  is  given  by 

ls(!s,z)=zG(2i^)^:f(Ks,z)  (C-16) 


where 


and 


20  =  -JT  lz  =  z  |Kg  -k| 


(C-18) 


The  dyadic  y  is  defined 


—  =  2lT  1  Cn^"  [(—  *  — )  (— *  ~)  "  ?  1  (C-19) 


To  a  good  approximation  y  can  be  written^ ^ 


2= 


(C-20) 


provided  waves  travel  nearly  normal  to  the  surface.  For  the  case  of  an 
ordinary  pump  wave,  an  ordinary  signal  wave,  and  an  extraordinary  sum- 
frequency  wave,  the  growing  wave  in  equation  C-16  becomes 


£sx(«s,  z)  =  2*i  ^  z  G(2  if )  P,  (Jt„,  z) 


X  '  s 


(C-21) 


Equations  C-16  and  C-21  can  be  written  in  more  familiar  fornr  3  by  noting 


G(2ilt)=e-i'i^l 


(C-22) 
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b.  EVALUATION  OF  MISMATCH  FACTOR 


Before  equation  C-16  can  be  applied  it  is  necessary  to  express 
the  mismatch  factor  0  in  terms  of  wave  vectors.  Equation  C-15  gives  that 

to 

0=  C1  (Ks  "  K)  '  N  (C-23) 

As  shown  in  Figure  C-2  the  wave  vector  K  of  free  waves  is  constrained  to 

to 

g 

lie  on  the  wave  vector  surface  -^-n  (s)  for  the  sum  frequency  to  ;  the  wave 

L  s 

vector  Kg  of  polarization  wave  is  given  by  Kg  =  K.  +  and  is  arbitrary. 

We  shall  denote  the  wave  vector  of  polarization  wave  in  the  matching  direc¬ 
tion  by  Km-  We  shall  assume,  for  simplicity,  that  the  surface  Normal  N  is 
a  matching  direction  and  also  that  is  in  the  direction  N.  Figure  C-2 
shows  geometrically  the  mismatch  relation  expressed  by  equation  C-23.  We 
have  taken  the  z  direction  to  be  along  N,  the  optic  axis  in  the  xz  plane,  and 
the  y  direction  normal  to  the  plane  of  figure.  The  tangent  plane  to  the  wave 
vector  surface  at  is  shown,  and  the  nor.nal  to  this  plane  makes  the 
angle  C  with  the  z  axis.  It  should  be  noted  that  the  angle  C  represents  aniso¬ 
tropy  of  the  crystal  at  the  frequency  us  .  The  anisotropy  angle  C  is  given  by: 

s 

—  tan  £  =  u?{9  )  cos  9  sin  9  (  2  ~  — (C-24) 

m  m  m\ne  nQc  } 

where  9^  =  angle  between  the  optic  axis  and  the  phase  match  directon 

n(0  )  =  index  of  refraction  along  at  u> 
m  -M  s 

ng  =  extraordinary-wave  index  of  refraction  at 

n  =  ordinary-wave  index  of  refraction  at  w 
o  s 
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This  angle  also  governs  the  double  refraction  effect,  in  which  both  the 
ordinary  and  extraordinary  waves  are  at  the  same  frequency.  In  this  analysis 
it  was  assumed  that  the  signal  field  and  the  pump  field  are  both  polarized 
as  ordinary  waves  and  the  sum-frequency  field  is  polarized  as  extraordinary 
waves.  The  mismatch  factor  0  according  to  equation  C-23  is  the  distance 
from  Kg  to  the  wave  vector  surface  measured  parallel  to  N.  This  is  apparent 
from  the  boundary  condition  that  requires  continuity  of  transverse  wave  vectors 
of  both  forced  and  free  waves  across  the  boundary  surface. 


tc  =  K  =  K  +  K 
s  -p  ~i 


(C-25) 


This  condition  assures  also  the  continuity  of  an  angular  spectra  with  spatial 
angular  frequency  across  the  crystal  boundary.  It  is  somewhat  tedious 
but  can  be  shown  that: 
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where 
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Since  the  crystal  is  usually  nearly  isotropic  and  the  quantity  in  the  square 
bracket  is  second-order  in  (K  -  J£)  and  smaller  than  the  linear  quantity, 
it  is  proper  to  use  simpler  form 


K  -  K  -  <  C 
sz  M  sx 


(C-27) 


We  have  used  tan  C  =  C  in  equation  C-27  since  C  <  <  1. 

c.  RESULTANT  SUM- FREQUENCY  FIELD 

The  sum-frequency  angular  spectrum  component  £  (H  ,  z)  given 

s  s 

by  equation  C-16  is  generated  by  a  plane  wave  polarization £(<  ,  C),  which 

s 

is  a  Fourier  component  of  total  polarization  wave  P(r).  The  resultant  sum- 

frequency  field  produced  by  P(r)  is  Fourier  transform  of  &  (K  ,  z) 

s  s 

Es(l)  =  J  6s(«s.z)  e  s  £dSg  (C-28) 

Substituting  equation  C-16  into  equation  C-28  the  resultant  field  can  be  written 

E(r)=z  f  dv  fl  :Jp(K,  z)  e'2l0T?  e"  '  -^dK  (C-29) 
s  jQ  j 

where  use  has  been  made  of  integral  representation  in  equation  C-17. 

o 

For  paraxial  assumption  the  quadradic  term  *S/2KM  can  be 
neglected.  It  then  follows  from  equation  C-29  and  inverse  Fourier  integral 
representation  of  £j£,  z)  that 
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-s(-  =^2  d  n  I  d  £<  e’-s '  s  -  KM  +  £  *sx> 

x/dp'y  :P(r')e'-s  -  (C-30) 

The  integration  over  k  gives  the  factor 

— s 

2 

(2tt)  6  (x  -  x'  -  t)  {  z)  5  (y  -  y') 

The  integration  over  p'  gives 

5g(£)  =  z  IQ  d  £  e  17?KM^  y  :  P  (x  -  *?£z,y,  z)  (C-31) 

Then  the  field  at  the  exit  surface  of  the  nonlinear  crystal  (z  =  D)  can  be 
written 

Eg(p,  D)  =  D  JQ  drje  17?KMD  y  ;  P(x-  t?£D,  y,  D)  (C-32) 

If  we  let  (1  -  rj)  D  =  z  and  make  z  as  a  dummy  variable  of  integration,  equa¬ 
tion  C-32  becomes 

Eg(p,  D)  =  e  iKMI')  /QD  da  eiKMZ  y  :  P(  x  -  £D  +  f  z,  y,  D)  (C-33) 

This  expression  shows  that  the  resultant  sum -frequency  field  is 
produced  by  constructive  interference  of  the  waves  radiated  from  each  dif¬ 
ferential  slab  of  polarization  in  the  material.  It  is  this  concept  that  was  used 
in  the  geometric  optics  theory.  Double  refraction  (or  anisotropy  effect)  of  the 
sum-frequency  electric  field  in  included  in  equation  C-33  through  the  param¬ 
eter  (x-CD  +  rz).  Now  refer  to  Figure  C-3  the  line  at  angle  C,  which  is  nor 
mal  to  the  tangential  plane  to  the  wave  vector  surface,  represents  the  direc¬ 
tion  of  energy  flow  of  the  sum-frequency  field(21).  It  is  intuitive  to  expect  that 
the  path  of  integration  to  obtain  resultant  field  should  be  along  the  energy 
flow  direction  rather  than  the  z  direction.  To  understand  the  implication  of 
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FIGURE  C-3. 

ILLUSTRATION  OF  DOUBLE  REFRACTION  OF  SUM 
FREQUENCY  WAVE 

the  parameter  x-^D  +  U,  let  us  denote  points  in  the  polarization  wave, 
the  source  points,  by  (x^  y  ,  z^,  and  points  in  the  sum-frequency  wave, 
the  observer  points,  by  (x2,  y2,  zg).  We  regard  equation  C-33  an  an  integral 
over  source  points  giving  the  sum -frequency  field  at  an  observer  point.  We 
shall  measure  the  source  points  from  the  coordinate  origin.  We  define  x^ 
for  the  observer  points  in  such  a  way  that  x  =  0,  y_  =  0  lies  on  the  line 
making  angle  £  with  the  z  axis.  Then  and  the  path  of  integration 

over  the  source  points  is  specified  by  x^  =  x^  +  £  z^,  y^  =  y2  =  y,  with  x^ 
hold  constant  in  the  integration  since  it  refers  to  the  observation  point.  Then 
when  we  set  z  =  D,  we  obtain  the  integration  parameter  x  =x-£D  +  £z. 

La  X 

This  can  be  interpreted  that  the  sum -frequency  field  at  the  points  (x,  y,D) 
is  obtained  by  integration  along  the  energy  flow  line. 

The  angular  spectrum  of  the  sum -frequency  field  at  the  exit 
surface  of  the  crystal  is  simply 


where 


•  aKD 

<5  (k  ,D)=De'1  2  y 

— s  — s  = 


*K=  K  -  K  -  £  k 
sz  n 


(C-34) 


C-35) 


Equation  C-34  indicates  that  the  angular  spectrum  at  the  exit  surface  is  not  a 
summation  over  the  contributions  of  the  Fourier  components  of  the  polariza¬ 
tion  in  the  differential  slab. 
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3.  FOURIER  SPACE  CONFIGURATION  (TYPE  III)  SYSTEM 

a.  POINT  SPREAD  FUNCTIONS  AND  OPTICAL  TRANSFER  FUNC¬ 
TION  OF  A  PLANE  WAVE  PUMPED  INFINITE  APERTURE 
FOURIER  CONFIGURATION  SYSTEM _ 

Diffraction  due  to  apertures  is  neglected  in  this  discussion  in 

order  to  separate  image  properties  of  the  upconversion  process  from  the 

diffraction.  Lens  apertures  (the  nonlinear  material  aperture)  and  the  pump 

beam  diameter  are,  therefore,  assumed  to  be  infinitely  large. 

A  simplified  schematic  of  a  Fourier  space  configuration  upcon- 

verter  system  is  shown  in  Figure  C-4.  Lens  L1  is  positioned  such  that  the 

object  and  the  front  surface  of  the  nonlinear  material  are  respectively  in 

the  front  and  back  focal  planes.  The  front  local  plane  of  lens  L0  coincide 

with  the  exit  surface  of  the  nonlinear  crystal.  Upconverted  real  image  then 

forms  in  the  back  focal  plane  of  lens  L 2»  Lenses  Lj  and  Lg  therefore, 

form  a  Fourier  transform  pair.  The  upconversion  takes  place  in  the 

Fourier  transformed  space  of  the  object  waves;  thus  the  name  Fourier  space 

configuration  system. 


For  the  plane  wave  pumped  upconverter,  it  can  be  shown  that 
the  angular  spectrum  of  upconverted  image  fields  in  air  beyond  the  non¬ 
linear  crystal  (z  >  +  D/2)  is  related  to  the  Fourier  component  of  object 


fields  over  the  front  surface  of  the  crystal  by 
£s  (*g>  z)  ~  £.  (*s,-D/2)  Sine 

e^sz^saz*2-0/^ 


■i  AKD/2 


(C-36) 
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FIGURE  C-4.  IMAGE  UPCONVERTER  ARRANGED  IN  FOURIER  SPACE 
CONFIGURATION 
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where 


Sine  x  =  (sin  x)/x 
D  :  material  thickness 

k  -Ik2  2 1  i/2 

Ksz^s  -  *  s  | 


z  component  of  wave  vector  at  X  in  the 
material  ® 


Ksaz  '  z  comP°nent  of  wave  vector  at  A„  in  air 


* 

1 

i 

I 

I 


A  K  (  Ksz  -  Kj^j  -  C  + 


sx 


2K 


)  :  phase  mismatch  factor 


M 


KM  :  wave  vector  at  along  the  phase  matched  direction 
C  :  anisotropy  angle 

The  angular  spectrum  £.  (  s  ,  -  D/2)  is  a  Fourier  transform  of  the  object 
field  distribution  over  the  plane  z  =  -D/2.  The  object  field  distribution  in 
the  plane  z  =  -D/2  is  a  Fourier  transform  of  the  object  field  distribution  over 
the  front  focal  plane  of  lens  since  the  plane  z=-D/2  coincides  with  the 


back  focal  plane  of  lens  L 


(14) 


1 


K. 


_i _ ra  j  '  \ 

E((p,  -D/2)  =  ^-|E.  (£.,  -2£l-D/2)e  ‘  dg  (c.37) 


and 


(C-38) 
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where 


f j  =  focal  length  of  lens  Lj 

_p^=  transverse  position  in  the  front  focal 
plane  of  lens 

__p  =  transverse  position  in  the  back  focal  plane 
(z  =  -D/2)  of  lens  Lj 

Substituting  Equation  (2)  into  Equation  (3)  and  carrying  out  the 
integration  we  obtain 

-  D/2)  -  E.  (-fj  /K.a,  -2f1  -  D/2)  (C-39) 

Constant  Amplitude  factor  and  fixed  phase  factors  have  been 
omitted  in  equation  C-39  since  they  do  not  affect  the  result  in  any  significant 
way.  Equation  C-39  indicates  that  the  angular  spectrum  of  field  distribution 
over  the  back  focal  plane  of  a  nondiffracting  thin  lens  has  the  same  functional 
form  as  the  field  distribution  over  the  front  focal  plane  of  the  lens.  This 
result  is  no  surprise  if  one  realizes  that  they  are  related  by  double  Fourier 
transform. 

The  field  distribution  of  an  upconverted  image  formed  in  the 
back  focal  plane  of  lens  L„  has  the  same  functional  form  as  £  {  k  ,  -D/2) 
since  the  plane  z  =  +D/2  and  the  upconverted  image  plane  are  two  focal 
planes  of  lens  Lg. 


C-21 


Es<fs’2f2  +  D/2>~  4(KsaPs/f2’D/2> 


(C-40) 


Combining  equations  C-36,  C-39,  and  C-40,  we  obtain  an  expression  for  the 
upconverted  image  field  distribution 


K  L 


Es  ^s’  2  f2  +  D//2^  ~  Ei  ("  K"“f  -£i»  "  2fl  "  D//2^ 

ia  2 


a  Kn  _i  -  K  D 
Sine  (-—^)  e  2  sz 


(C-41) 


where  mismatch  factor  AK  and  the  z  component  wave  vector  K  in  equa- 

SZ 

tion  C-41  are  evaluated  by  replacing  jc  with  K  p  / f  and  are  given  by 

S  S3.“S  ct 


AK  =  Ks-KM-«KsaVf2  +  Tf2-Ps2  <44>  <C-42> 

2  f  g  B  M  s 


Ksz2=Ks2-<Ksa'’/f2>2 


(C-43) 


If  crystal  anisotropy  is  neglected  the  mismatch  factor  AK  is 


AK  =  - 


The  function  sine 


f2 


K  D  K 


(C-44) 


4KMKi 


indicates  that  the  amplitude 


of  an  upconverted  image  field  of  a  uniformly  illuminated  object  field  will  not 
be  uniform  but  will  decrease  with  transverse  distance  pg.  It  imposes  a 
limitation  on  the  size  of  the  upconverted  image  but  will  not  limit  the  image 
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resolution.  If  we  define  the  size  of  an  upconverted  image  as  equal  to  the 
transverse  distance  pg,  where  the  first  zero  of  a  sine  x  function  occurs, 
the  radius  of  an  upconverted  image  is  expressed  by 


a  =  2— 
s  K 


sa 


n  K™K- 
M  1 


D  K 


(C-45) 


Equation  C-41  also  shows  that  the  transverse  magnification  of  a  Fourier  space 
configuration  upconverter  system  is 


M  =  -  Kia  *2  -  Xg  *2 

‘  Ksafl  '  Vl 


(C-46) 


Equation  C-46  can  also  be  obtained  by  the  geometric  optics  analysis.  The 
limiting  size  of  an  object  that  can  be  upconverted  can  be  easily  found  by 
combining  equations  C-45  and  C-46. 


vVh  =2 


'rr  K.  _  K. 
M  1 


K. 


la 


DK 


(C-47) 


Noting  that  the  object  size  a^^  is  related  to  the  acceptance  angle  0ia(FOV)  of 
the  upconverter  by 


2ai  ■  fl 


e . 

ia 


we  obtain 


2  X  X  .  n .  n 
p  l  l  sM 


\  1/2 


X  n  D 
s  p 


(C-48) 


0 


li 


a 

i. 


I. 

T 

«  » 


For  a  noncriticaliy  phased  matched  proustite  unconversion  from  10. 6  urn  to 

i  \1/2  . 

D  rad' 


0. 967  um,  we  have  0.  =0.157 
^  ’  la 


To  illustrate  the  points  presented  in  the  previous  discussion  let  us  consider 
a  case  where  the  object  field  distribution  is  a  sine  function  as  shown  in 
Figure  C-5A.  The  field  distribution  of  upcon verted  image  formed  in  the  back 
focal  plane  of  lens  L 2  is  shown  in  Figure  C-5B.  The  period  of  the  sine  function 
image  field  is  magnified  by  a  factor  and  its  amplitude  envelope  follows 
the  sine  function.  The  image  is  inverted  as  shown  in  Figure  C-5B  so  that  an 
object  point  -y^  is  imaged  at  yg^  =  -M^  y^.  The  resolving  power  is  un¬ 
limited  since  the  image  field  period  is  linearly  proportional  to  the  object 
field  period.  The  object  size  that  can  be  upconverted  is,  however,  limited. 

The  size  limitation  is  caused  by  the  phase  match  requirement  rather  than 
the  crystal  aperture  limitation. 

The  point  spread  function,  which  is  a  response  to  a  point  object, 
is  easily  obtained  from  equation  C-41.  The  field  distribution  of  a  point 
object  at  (x.q,  y.  )  is  expressed 

Ei  <£i>  ~  6  <xi  -  xio)  •  6!yi  -  yio> 


Then  from  equation  C-41,  the  point  spread  function  becomes 
Hs  <?s>  ~  6  <xs/Mt  ■  xio>  6<ys/Mt  •  yio> 


-i  (AK/2  -  Ksz)  D 


Sine 

fAKE^ 

2 

/ 

-  ~  Ksa  Mt  Pio  /f2 


(C-49) 
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B,  UPCONVERTED  IMAGE  FIELD  DISTRIBUTION 


FIGURE  C-5.  FIELD  DISTRIBUTION  OF  OBJECT  UPCONVERTED  IMAGE 
FOR  THE  FOURIER  SPACE  CONFIGURATION 


where 


k  =  Ksa£s  indicates  that  the  quantity  inside  the 
bracket  is  evaluated  at^c  =  K  p  /f  as  shown  in  equations  C-42  and  C-43. 

S3.  —  S  ct 

The  coherent  transfer  function,  which  is  the  Fourier  transform 
of  the  point  spread  function,  becomes 


2  -i  M  £  •  p 

hs  <*s>  ~  Mt  e 


10 


Sine 


A  KD 


i 

[a,k-k  1 

D 

[  2  szi 

(C-50) 


=  K  M,  p.  /f. 
sa  t  _io  2 


Amplitude  of  h  (  k  )  depends  on  the  object  location  and  thus  object  size 


limitation  is  implied  even  for  the  infinite  aperture  Fourier  configuration 

14 

vpconverter.  The  optical  transfer  function,  which  is  defined  as, 


ws->  ■ 


JJhs 


(C-51) 


d  k  ' 


becomes 


Hs  (^)=  e 


i  k  •  p 


10 


(C-52) 


Equation  C-52  again  indicates  the  property  of  unlimited  resolving  power  since 
the  amplitude  of  optical  transfer  function  (Modulation  transfer  function)  is 
independent  of  angular  frequency  k  . 


b.  POINT  SPREAD  FUNCTION  OF  A  GAUSSIAN  BEAM  PUMPED 

FOURIER  SPACE  CONFIGURATION  SYSTEM _ _ 

The  effect  of  pump  beam  amplitude  distribution  on  the  imaging 
characteristics  of  the  Fourier  space  configuration  image  upconverter  is  con¬ 
sidered.  Diffraction  due  to  the  lens  apertures  and  the  crystal  birefringence 
are  neglected.  It  is  also  assumed  that  the  physical  size  of  the  nonlinear 
material  is  large  compared  to  the  pump  beam  size.  The  optical  arrange¬ 
ment  is  shown  in  Figure  C-4. 

The  pump  beam  is  assumed  to  be  in  the  TEM  mode  The 

oo 

Gaussian  window  at  the  entrance  surface  of  the  nonlinear  crystal  is 
expressed  by* 


2/2 

Ep(o  ,  -D/2)  ~  e  /a 


The  angular  spectrum  of  pump  field  distribution  is 


-D/2) 


(C-53) 


(C-54) 


The  point  object  in  the  front  focal  plane  of  the  object  lens  is  expressed  by 


E  (o  ,  -2f  -D/2)  ~  6  (o.  -  P.  ) 
1—1  1  -1—10 


(C-55) 


*Constant  amplitude  and  phase  factors 
analysis. 


are  again  omitted  in  the  following 
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The  angular  spectrum  in  the  IR  object  fields  in  the  entrance  plane  of 
the  crystal  is,  using  equation  C-39 

6j(  ^,-0/2)  -  E.MjS./K^,  -2^0/2)  -  6(  S,  +  K_a  p  ^/fj)  (C.B6) 


The  sum-frequency  polarization  P  is  obtained  by  taking  a  convolution  of 

s 

£.(*.,  -D/2)  and  £  (  k  ,  -D/2), 

l  '  ”i *  '  ’  P  P 

P  (  k  .  -D/2)  /  &.  (  k.  ,  -D/2)  £(*-*,  -D/2)  d  < . 

S  -s’  '  '  j  1  l'  p  s  p’  '  '  1 


-a2  (  k  +K.  p.  /L)2/4 
s  ia-10  1  / 


(C-57) 


The  angular  spectrum  of  upconverted  image  fields  at  the  exit  surface  of  the 
crystal  is 


£  (  « 
s  s 


+D/2) 


•a2  (  «  +K. 

'  s  la 


VV2/« 


Sine 


‘  \  -iD(K  +^) 

AKD  V  sz  2  J 

.  2  J 

(C-58) 


where  AK  and  K  are  mismatch  factors  and  the  z  component  of  wave 
sz 

vector  in  the  material  for  the  angular  spectrum  with  the  angular 
frequency  k  . 

Field  distribution  of  the  upconverted  image  in  the  back  focal  plane  of  the 
imaging  lens  is,  by  equation  C-40 

Es(-°s'2[2+D/2)~  es(K6a-°s/f2’  +D/2' 


C-28 


0  /  K 
2  sa 


sa  la 

7 —  +  "7 —  o. 

f2  S  fl  _1° 


fsinc  e-‘D<Ksz*iK/2) 

L  J  *  = 


K 

5  =-S^  p 
s  fo  S 


(C-59) 


The  symbol 


v  indicates  again  that  AKand  K  are 
■^sa  sz 


—  c  u 
S  I9  -S 


evaluated  by  replacing  it  with  K  o  /f0  as  shown  in  equations  C-42  and  C-43. 

o  Sa  S  A 

Equation  C-59  is  the  point  spread  function  of  Gaussian  beam 
pumped  Fourier  space  configuration  image  upconverter.  The  upconverted 
image  is  no  longer  a  point.  The  sine  function  again  limits  the  FOV  of  the 
system  and  describes  the  amplitude  envelope  of  the  image  plane  in  the 
similar  way  as  the  plane  wave  pumped  case.  The  Gaussian  amplitude  factor 


aK. 
_ la 

2L 


fi  K 

1  sa 

r~v -  0  +  o. 

f2  Kia  "s  -,0 


indicates  that  transverse  magnification  factor  and  the  image  spot  size  are 


given  by 


Mt  =-°s/-i  =  -f2Kia/flKsa 


(C-60) 


W  =  2  f0/a  K 
s  2/  sa 


(C-61) 


Equcation  C-60  is  identical  to  equation  C-46.  It  implies  that  the  Gaussian 
amplitude  distribution  does  not  introduce  distortion.  Combining  equations 
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C-60  and  C-61  and  dividing  the  object  spot  size  by  the  focal  length  of  the 

object  lens  we  obtain  expressions  for  the  resolvable  IR  object  spot  size  and 
angle; 


wi  =  Vi/™ 

(C-62) 

y.  =  X.  /  va 

1  i/ 

(C-63) 

Since  Rayleigh  resolution  criteria  cannot  be  applied  to  the  present  case, 
the  1/e  amplitude  is  arbitrarily  chosen  as  the  resolution  limit. 

From  above  discussion  one  may  conclude: 

•  POV  and  the  intensity  envelope  of  image  plane  may  not 
be  affected  by  the  pump  beam  intensity  distribution. 

•  The  optical  transfer  function  of  the  Fourier  space 
configuration  system  depends  heavily  on  the  pump  beam 
intensity  distribution.  A  multimode  pump  beam  (diverging 
pump  beam),  therefore,  degrades  resolution  of  the  Fourier 
space  configuration  upconverter. 

4*  MAGE  UPCONVERFfQN  IN  IMAGE  SPAPT? 

Image  upconversion  with  Type  II  optical  system  ( Figure  C-6) 
is  considered  first  before  analyzing  Type  I  system.  As  shown  in  Figure  C-6. 
the  Type  II  system  does  not  use  an  IR  objective  lens.  For  a  distant  IR  object, 
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this  optical  system  is  very  similar  to  the  Type  III  system.  Diffraction 
due  to  the  crystal  and  other  apertures  and  material  birefringence  are  ignored 
in  the  following  discussion. 

The  pump  beam  intensity  distribution  at  the  crystal  surface  is 
assumed  to  be  Gaussian 

2  2 

Ep  (p,  -D/2)  =  Ep  e"p  /wp  (C-64) 

Its  angular  spectrum  is  expressed  by 

E  w  2  -w  2  «  2/4 

gp(«p,  -D/2)  =  ■  e  p  p  (C-65) 


If  the  angular  spectrum  of  IR  object  field  on  the  front  surface  of  crystal  is 
denoted  by  (K.,  -D/2),  the  angular  spectrum  of  the  sum-frequency  field 
in  the  space  z  >  D/2  is  given,  by  equations  C-21  and  C-22, 


6  (K, «>  eiD(Ksz  -  ^  Si»cUf& 

S  ' 


(C-66) 


exp  [i  Kgaz  (z  -  D/2)]  J  djc  #  6.  (I#  -  D/2)  exp  [-  w  2  (j£g  -  K#)2/4] 


where: 


ts,  tp,  t.  =  transmittance  at  the  crystal  surface  for  Xg,  Xp,  and  X^, 
respectively 

Sine  X  =  (Sin  X)/X 
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D  =  material  thickness 


II 


i  1 


0 


K  -  z  component  of  wave  vector  at  X  in  the  material 

sz  s 

K  =  z  component  of  wave  vector  at  X  in  air 

odZ  g 


K 


AK  =  K  -  K  -CK  + 


sz  m  sx  2K. .  I 

M 


:  phase  mismatch  factor 


KM  ”  wave  vector  at  Xc  in  the  material  along  the  phase- match 


direction 


il 


C  =  anisotropy  angle 
y  =  nonlinear  susceptibility 


If  a  point  IR  object  located  at  a  point  ( p^,  -s)  is  expressed  by 


E 


E.  (p.,  -s)  =  E.  6  (p.  -  p.  ) 
1  -1  1-1  -10 


(C-67) 


its  angular  spectrum  at  z  =  -D/2  (that  is,  on  the  crystal  front  surface)  is 
expressed  by 

E 

*1  (Si.  ■  E>/2)  - - 1-W  exp  i  [K  (s  -D/2)  -  (K  •  p  +  0jt  2)]  (C-68) 
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K.  =  wave  vector  at  X.  in  air 
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a  =(s  -  D/2)/2  K. 
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Substituting  equation  C-68  into  equation  C-66  and  carrying  out  integration, 
we  obtain  for  z  >  D/2 


„  co  E.E  w  t  t.t  vD  . 

(£„>  z)  =  i  —  1  P  P~2  P  1  S - sinc  (— o^Jexp  [i  0  (if,  z)] 

2tt  Cn  (w  +  i4  a )  S 

s  '  p  ' 


s  •  s 


(C-69) 
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where  the  phase  factor  is  expressed  by, 


I <IS.  z)  =  [Ksa  (z  -  D/2)  +  Kia  (s  -  D/2)  +  D  (K  +  K  )/2]  +  ^ 


sa 
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z  -  D/2  ,  s  -  D/2  D  ,  1  L 

v-tr  +  v> 
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K.  '  2  VK._  K' 
sa  la  Ms, 


2  sx 
(C-70) 


The  quantity  within  the  first  bracket  of  equation  C-70  is  independent  of  k 

s 

and  represents  a  constant  phase  shift  for  all  angular  spectra.  The  factor 
I  j*d  \ 

\  "2~  *sx  I  arises  from  crystal  birefringence.  Its  effect  on  the  imaging 
properties  will  be  discussed  in  detail  later,  but  it  will  be  ignored  for  the 
present  discussion.  The  stationary  phase  condition  that  is  required  for  the 
reproduction  of  clear  image,  the  phase  factor,  requires  thatjJ  should  be 
independent  of  spatial  angular  frequency.  This  condition  is  used  to  determine 


the  relationship  between  two  conjugate  planes.  Equating  the  coefficient 
o 

of  k  to  zero,  we  can  solve  for  the  image  plane  distance 
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The  angular  spectrum  of  virtual  upconverted  image  field  is  now  given  ignoring 
a  constant  phase  factor  by 
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The  exponential  decaying  factor  in  equation  C-72  is  resulted  from  the 
Gaussian  window  due  to  the  pump  beam  intensity  distribution.  The  exponen¬ 
tial  decaying  function  becomes  unity  if  the  pump  beam  is  infinitely  large 

(w  =  a.)  that  is,  plane  wave  pump.  The  exponential  phase  function  arises 
P 

from  the  off-axis  object  location. 

The  point  spread  function  of  the  virtual  upconverted  image 


is  obtained,  for  negligible  birefringence,  by  using  a  relation 
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and  taking  a  Fourier  transform  of  equation  C-72 
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The  point  spread  function  of  Type  I  system  can  easily  be  obtained 
from  equation  C-73  by  equating  a  to  zero.  This  can  be  understood  if  one 
realizes  the  geometry  of  the  Type  I  system.  In  the  Type  I  system,  the 
objective  telescope  is  used  to  obtain  a  real  IR  image  of  the  object  in  the  plain 
containing  the  pump  window.  The  factor  a  is  equal  to  zero  since  s  =  D/2 
for  Type  I  arrangement. 

The  point  spread  function  of  Type  I  system  is,  therefore,  after 
the  integration, 
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The  exponential  decay  factor  is  a  Gaussian  and  describes 
the  envelope  of  upconverted  image  field  amplitude  distribution.  The  sine- 
integral  function  describes  the  image-object  relations.  It  is  obvious  from 
elation  C-74  that  the  upconversion  transverse  magnification  is  unity  for 
Tjpe  I  system.  Taking  the  first  root  of  sine- integral  function  as  the  limiting 
resolution,  we  obtain  an  expression  for  the  linear  resolution  for  Type  I  system. 


5.  PARAMETRIC  IMAGE  UPCONVERSION  WITH  A  MULTIMODE 
LASER  PUMP 

Parametric  image  upconversion  pumped  by  plane  waves  and 
Gaussian  beams  have  been  considered  in  proceeding  sections.  It  was  shown 
that  diffraction-limited  resolution  can  be  obtained  under  certain  conditions. 

To  make  full  use  of  the  available  pump  laser  power,  it  would  be  necessary 
to  operate  a  pump  laser  in  multimode  pattern.  The  imaging  properties  of 
a  Type  II  system  using  a  multimode  pump  laser  are  considered  first  in  this 
section. 

The  spatial  frequency  response  (angular  spectrum)  of  the  up- 
converted  image  field  is  considered  for  two  cases.  In  the  first  upconverter 
system,  the  nonlinear  material  is  pumped  by  the  near  field  of  a  TEM^^ 
mode  laser  beam.  The  crystal  in  the  second  system  is  placed  in  the  near 
field  of  a  toroidal  mode  beam.  The  shape  of  the  spatial  frequency  response 
resembles,  as  expected,  that  of  pump  field  distribution.  The  result  suggests 
possible  spatial  filtering  technique  for  the  upconversion  system.  Quantitative 
evaluation  of  the  field  distribution  of  upcon verted  image  was  the  attempted 
at  this  time  since  it  requires  an  extensive  electronic  computer  effort. 

Imaging  properties  of  the  Type  I  and  III  systems  pumped  with  a  TEMj(1 
mode  and  a  troidal  mode  pump  beam  are  then  considered  qualitatively. 


Consider  first  the  case  of  a  TEM^j  node  pump  laser  beam. 
Assume  that  the  material  is  placed  in  the  near  field  (the  waist)  of  the  beam. 
The  field  distribution  of  the  pump  beam  can  be  considered  as  an  optical  win¬ 
dow  for  IR  object  waves.  The  near-field  distribution  of  a  TEMQ1  mode  beam 
is  expressed,  in  a  cylindrical  coordinate  system  as, 
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where 


p  =  radial  coordinate 
a  =  beam  radius  of  fundamental  mode 

Unity  amplitude  is  assumed  in  equation  C-76  for  simplicity.  The  ang.  <ar 
spectrum  of  a  TEM^  field  distribution  is  the  Hankcl  transform  of  field 
distribution  g'ven  by  equation  C-76, 
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It  can  be  shown  that  6  (<),  neglecting  a  constant  multiplicative  factor, 

becomes 
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Neglecting  double  refraction  effects  and  following  the  procedure  developed 
previously,  it  can  be  shown  that  the  equation  relating  the  image  position 
and  the  object  position  is  identical  to  that  of  Gaussian  beam  pumped  case. 
The  angular  spectrum  of  upconverted  field  distribution  over  the  image  plane 
for  an  on-axis  point  object,  again  neglecting  constant  multiplicative  factors, 
is  given  by: 
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where 

a  =  (s  -  D/2)/2v 


The  (sin  x/x)  function  arises  again  from  phase  mismatch  for  a  given  angular 
spectrum.  The  attenuation  factor  expressed  by  the  integration  in  equation  C-79 
is  due  to  the  transverse  mode  structure  of  the  pump  laser  beam.  The  attenua¬ 
tion  factor, 
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can  h(?  rewritten  as  Hankel  transform. 
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Evaluating  the  integration,  we  obtain 
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If  the  condition  a  »  2a  is  satisfied,  equation  80  becomes: 
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Equations  C-80  and  C-81  indicate  that  the  image  spectrum  at  high  spatial  fre¬ 
quencies  of  the  TEMqj  beam  pumped  image  upconvcrtcr  system  is  attenuated 
more  than  that  of  Gaussian  beam  pumped  system. 

We  consider,  next,  an  upconvcrtcr  system  pumped  by  a  toroidal 
mode  laser  beam.  The  toroidal  mode,  the  simplest  multimode  distribution, 
considered  here  is  pr?*iuccd  by  subtracting  equal  amplitude  TEM(|1  mode  field 
from  TEMqo  mode  field.  The  near-field  distribution  of  this  mode  is  expressed 
by: 

Ep(p)  B  EOOW  "  E01  W  ■  2  exp  )  (C-82) 
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The  angular  spectrum  of  the  toroidal  field  distribution  is: 
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It  can  be  shown  that  the  mode  attenuation  factor  for  t  ic  angular  spectrum  of 


the  upconvcrted  image  field  produced  by  this  system  is 
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Again  for  the  condition  a  >>  2a,  the  attenuation  factor  approaches 
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Equations  C-85  and  C-8G  indicate  that  a  relative  enhancement  at  high  spatial 


frequency  takes  place  with  the  toroidal  mode  pumping.  The  dc  frequency 


spatial  component  Is  decreased  and  the  peak  frequency  response  shifted  to 


higher  frequencies. 


The  above  results  strongly  suggest  that  spatial  filtering  can  be  per¬ 


formed  by  properly  shaping  the  pump  beam.  Many  techniques  developed  for 
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optical  data  processing  and  holography  can  be  applied  directly  to  the  para¬ 
metric  image  upcooverter  if  the  pump  beam  shape  is  predictable  and  stable. 

The  above  computational  technique  can  be  applied  to  any  TEM^ 

(n  =  0,1,2,...)  mode  pump  beam.  It  is,  however,  expected  that  the  mode 
attenuation  factor  for  circularly  asymmetrical  modes  will  not  be  expressible 
in  a  closed  form. 

For  a  TEMQ1  mode  beam  pumped  Type  1  system,  the  attenuation 
factor,  expressed  by  equation  C-80  of  angular  spectrum  of  the  upconverted 
image  of  an  on-axis  point  object  becomes  unity  since  a  is  identically  zero. 

This  implies  that  the  factor  sine  (A  K  D/2)  will  be  the  dominant  factor 
affecting  the  resolution  of  the  upconverted  image.  The  resolving  power  of 
Type  I  system  is  therefore  independent  of  the  pump  beam  mode  structure. 

For  a  toroidal  inode  pumped  Typo  i  system,  the  attenuation  factor  equation  C-84 
becomes  identically  equal  to  zero  for  an  on-axis  point  object.  It  implies  that 
on-axis  objects  would  not  be  upconverted.  This  is  no  surprise  if  one  realizes 
that  the  on-axis  pump  power  is  zero  and  object  waves  interact  with  local  pump 
beam  in  the  Type  I  system.  It  is  obvious  from  above  observation  that  the 
envelope  of  upconverted  image  field  of  Type  I  system  would  have  the  same 
functional  form  as  that  of  pump  beam. 

It  is  clear  from  equations  C-80  and  C-86  that  the  angular  spectrum 

of  upconverted  field  for  Type  III  system  depends  heavily  on  the  mode  structure 

2 

of  pump  beam,  since  a  >  a  for  Type  III  system.  The  resolution  and  the 
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point  spread  function  of  Type  III  system  are,  therefore,  influenced  strongly 
by  the  pump  beam  mode  structure.  It  can  be  shown  easily  that  the  point 
spread  function  of  Type  III  system  pumped  with  a  uniform  truncated  plane 
wave  (a  pill-box  pump)  follows  a  J^(X)/X  function.  It  has  also  been  shown  in 
paragraph  3,  that  the  point  spread  function  of  TEMqo  mode  beam  pumped 
Type  III  system  is  a  Gaussian  function.  These  examples  suggest  that  the 
image  resolution  of  Type  in  system  will  be  degraded  by  the  pump  beam 
divergence. 

6.  EFFECTS  OF  CRYSTAL  BIREFRINGENCE 

In  this  paragraph  we  will  consider  effects  of  the  crystal  birefrin¬ 
gence  upon  imaging  properties  of  the  upconverter.  Negligible  aperture 
diffraction  and  plane  wave  pumn  are  assumed  in  order  to  isolate  the  birefrin- 
gent  effect  from  others. 

a.  IMAGE  SPACE  CONFIGURATIONS 

Omitting  constant  amplitude  and  phase  factors,  the  angular 
spectrum  of  upconverted  field  emanating  out  of  the  crystal  back  surface  i.s 
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where 
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The  electric  field  distribution  in  the  upcon verted  image  field  becomes 
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Equation  C-89  indicates  that  the  optimum  focusing  condition  is 
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is  identical  to  the  case  of  negligible  birefringence  system  and  defines  the 
image  plane. 

The  line  spread  function  of  a  line  object  parallel  to  the  x  axis 
will  be  considered  first.  The  object  field  distribution  and  its  angular 
spectrum  are 
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The  line  spread  function  is  then 
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Evaluating  the  integration,  we  obtain 
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where  Fresnel  functions  S(X)  and  C(y)  are  defined 
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The  expression  C-92  applies  equally  to  the  system  with  negligible 
birefringence.  This  result  indicates  that  the  image  quality  is  not  degraded 
along  the  y  direction  which  is  normal  to  the  plane  containing  optic  axis  C 
and  the  z  axis  (Figure  C-l).  The  line  spread  function  is  only  a  function  of 
normal  distance  (y„  -  y.  ). 
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The  line  spread  function  of  a  line  object  parallel  to  the  y  axis  is 
considered  next.  The  object  field  distribution  and  its  angular  spectrum  are 
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Rewriting  above  equation, 
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where 


Equation  C-94  indicates  that  crystal  birefringence  affects  the  upconverted 
image  in  following  ways : 

•  Introducing  change  in  the  field  distribution  from  that  of 
non-birefringent  case 
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•  The  field  distribution  is  a  function  of  (x  -  x.  -  C  D/2) 

s  10 

rather  than  x„  -x.  .  This  implies  that  the  x  coordinate 
of  upconverted  image  is  translated  toward  positive  x 
direction  by  an  amount  CD/2,  (Figure  C-7).  This  linear 
translation  becomes  obvious  if  one  realizes  that  energy 
flow  inside  the  material  is  along  the  anisotropy  direction 
rather  than  along  the  z  direction. 

Line  spread  functions  along  two  orthogonal  directions  expressed  by  equation 
C-92  and  equation  C-94  are  evaluated  for  a  1.06-fim  pumped  1  cm  long 
proustite  crystal  upconverting  10.6-um  objects  into  0.967-jim  images. 

The  results  are  plotted  in  Figure  C-8.  The  amplitude  of  the  spread  functions 
is  normalized  to  the  peak  amplitude  of  the  line  spread  function  parallel  to 
the  x  axis. 

Resolvable  image  line  width  along  the  y  direction  can  be 
determined  by  Rayleigh  criterion  since  the  sidelobes  of  the  spread  function 
parallel  to  the  x  axis  are  considerably  lower  compared  to  the  mainlobe. 

The  spread  function  parallel  to  the  y  axis  contains  several  nearly 
equal  amplitude  lobes  as  seen  in  Figure  C-8.  The  amplitude  of  these  lobes  is 
considerably  lower  than  that  of  the  mainlobe  of  the  x  direction  line  spread 
function.  The  width  of  these  lobes  is  also  narrower  than  that  of  the  x  direction 
line  spread  function.  Rayleigh  resolution  criterion  cannot  be  used  to  deter¬ 
mine  resolution  along  the  x  direction,  because  the  adjacent  lobes  are  nearly 
equal  in  amplitude. 


Resolution  degradation  along  the  x  axis  predicted  from  Figure  C-8 
has  not  been  observed  in  the  experiments  described  in  Appendix  A.  One  of 
possible  causes  for  not  observing  resolution  degradation  is  due  to  the  masking 
of  birefringent  effect  by  relatively  larger  diffraction  line  width  of  the  object 
lens. 

More  careful  measurements  of  line  spread  function  parallel  to 
the  y  axis  are  necessary  in  order  to  observe  the  effect  of  crystal 
birefringence. 


b.  FOURIER  SPACE  CONFIGURATION 


It  has  been  shown  in  paragraph  A  that  the  field  distribution  of 
upconverted  image  formed  in  the  back  focal  plane  of  lens  L 2  (Figure  C-4)  is 
related  to  the  object  field  distribution  over  the  front  focal  plane  of  lens 
by  repeating  equations  C-41,  C-42,  and  C-43 
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where  AK  and  K  are  given  by  equations  C-42  and  C-43 
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From  previous  equations  we  can  conclude: 


Magnification  along  two  orthogonal  transverse  directions 
is  not  affected  by  the  crystal  birefringence.  This 
implies  no  distortion. 

Different  amplitude  envelope  along  two  orthogonal  trans¬ 
verse  axes.  This  implies  different  acceptance  angle 
for  two  directions.  Acceptance  angle  in  the  yz  plane 
is  given  by  equation  C-48.  In  the  xz  plane  it  becomes: 
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No  degradation  in  resolving  power,  since  the  factor 
Sine  1  in  equation  C-41  only  affects  the  envelope  of 
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MATERIAL  DAMAGE  AND  LIGHT  ABSORPTION 
OF  PROUSTITE  AT  1. 06  pm 

1.  MATERIAL  DAMAGE 

The  question  of  material  damage  is  central  to  the  design  and 
operation  of  a  parametric  upconversion  system;  the  reason  being  the  large 
pomp  radiation  power  densities  to  which  the  nonlinear  materials  are  sub¬ 
jected  for  attaining  high  conversion  efficiency. 

Damage  thresholds  of  the  proustite  crystal  under  different  con- 
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ditions  have  been  reported  by  many  authors.  Fountain  et  al  reported  on 

damage  thresholds  of  proustite  at  1. 06  pm  subjected  both  the  CW  irradiation 

and  a  Q-switched  laser  pulse.  They  ascribed  damage  in  proustite  to  local 
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heating  caused  by  absorption.  Hanna  et  al  also  published  their  results  of 

Q-switched  laser  damage  of  proustite  crystals  at  1.06  pm,  10.6  pm,  and 

694  nm.  The  damage  thresholds  of  Q-switched  laser  irradiation  at  1.06  pm 

reported  by  Fountain  and  Hanna  are  comparable.  Hanna,  however,  excluded 
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surface  melting  due  ta  absorption  as  a  probable  damage  mechanism.  Lucy 

has  also  reported  damage  threshold  of  proustite  with  a  pulsed  Ruby  laser. 

2 

He  observed  significant  damage  at  peak  power  densities  of  10  kW/cm  with 
1-ms  pulse.  This  level  is  three  orders  of  magnitude  below  those  reported 
by  Hanna  with  a  Q-switched  Ruby  laser.  Accumulation  of  damage  is  apparent 
from  these  data. 
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The  three  cited  references  have  not  reached  a  clear  correlation 
between  damage  thresholds  and  bulk-absorption  coefficients.  They  disagreed 
on  a  probable  damage  mechanism. 

At  AIL  the  upconversion  system  has  been  operated  by  pumping 

the  proustite  crystal  with  CW  Nd:YAG  laser  radiation.  An  accumulation  of 

white  powder  on  the  exit  surface  of  proustite  crystal  has  been  noticed.  This 

23  25 

may  be  attributed  to  radiation  enhanced  surface  oxidation  ’  .  The  area 

and  shape  of  powder  accumulation  matches  the  1.06-pm  beam  cross  section. 

The  powder  accumulation  is  noticeable  after  a  few  tens  of  minutes  of  con- 

tinuous  irradiation  at  approximately  1  W/cm  .  Warner  indicated  the  same 

phenomenon  with  a  continuous  694-nm  irradiation  of  proustite.  He  reported 

that  the  surface  oxidation  becomes  visible  after  several  minutes  of  694-nm 

2 

irradiation  at  approximately  1  W/cm  . 

In  an  attempt  to  obtain  more  information  on  the  surface  oxidation 
of  proustite,  an  experiment  was  performed  by  placing  a  proustite  in  an  evac¬ 
uated  vessel.  The  vessel  was  evacuated  to  approximately  one  Torr  pressure 
and  the  crystal  was  irradiated  with  a  chopped  1. 06- pm  beam  of  80-percent 

duty  factor.  The  average  power  density  of  1. 06 -pm  beam  was  approximately 

2  2 
3  W/cm  with  a  temporal  peak  power  density  of  37.  5  W/cm  .  Surface  oxida¬ 
tion  was  not  observed  even  after  1/2  hour  of  irradiation.  This  level  of  ex¬ 
posure  is  considerably  higher  than  the  level  proustite  would  tolerate  in  the 
atmosphere.  The  result  indicates  that  Q-switched  laser  damage  threshold 
of  proustite  might  be  improved  if  the  proustite  is  operated  in  a  vacuum  or  in 
an  inert  gas  atmosphere. 
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2. 


LIGHT  ABSORPTION  OF  PROUSTITE  AT  1.  06  pm 


The  correlation  between  damage  thresholds  and  bulk  light 
absorption  is  not  clearly  understood.  There  are,  however,  data  that  in¬ 
dicate  that  the  damage  thresholds  and  surface  oxidation  levels  of  proustite 
are  lower  at  694  nm  than  at  1.  06  pm.  Lower  damage  threshold  at  694  nm 
may  be  attributed  to  higher  light  absorption.  Therefore,  it  has  been  found 
to  be  worthwhile  to  obtain  the  light  absorption  coefficient  of  proustite  at 
1.  06  pm,  since  bulk  absorption  of  proustite  at  1. 06  pm  has  never  been  re¬ 
ported  in  literature. 

A  calorimetric  measurement  was  done  of  the  absorption  coef¬ 
ficients  of  a  proustite  crystal  prepared  for  upcooversion  of  10. 6-pm  in¬ 
frared  radiation  to  967  nm.  *  The  crystal  is  oriented  such  that  the  optic 
axis  (C-axis)  makes  a  20-degree  angle  with  the  normal  of  the  (polished) 
front  surface.  The  crystallographic  X-axis  is  parallel  to  the  front  surface. 
A  linearly  polarized  chopped  1. 06-pm  beam  from  a  CW  TEMqo  mode 
Nd:YAG  laser  was  incident  normally  on  the  crystal  front  surface.  The 

average  power  densities  at  the  center  of  1. 03-pm  beam  were  approximately 
o 

5  W/cm  .  The  duty  factor  of  light  chopper  is  8.  5  percent  and  the  repetition 
period  is  short  enough  to  neglect  the  thermal  time  constant  of  the  material. 


*  Measurements  were  performed  with  an  absorption  meter  at  Holobeam  Inc. 
Paramus,  N.  J. 
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The  average  power  level  incident  on  the  crystal,  the  average 
power  transmitted  through  the  material,  and  the  average  power  absorbed  by 
the  crystal  were  measured.  The  proustite  was  mounted  on  a  Peltier  cooler. 
Electric  current  through  the  Peltier  cooler  was  adjusted  until  an  equilibrium 
was  reached  between  the  heat  generated  by  the  crystal  light  absorption  and 
the  cooling  rate  of  the  Peltier  cooler.  The  amount  of  light  power  absorbed 
by  a  material  is,  thus,  measured  by  observing  the  cooling  current  of  the 
calibrated  Peltier  cooler. 

Absorptions  of  the  ordinary  wave,  which  is  polarized  normal 
to  the  crystal  C-axis,  and  of  the  extraordinary  wave,  which  is  polarized 
parallel  to  the  C-axis,  were  measured.  Reflection  losses  at  the  surfaces 
were  calculated  from  the  Fresnel  reflection  formula  sin~e  the  index  of  re¬ 
fraction  of  proustite  is  accurately  known.  Table  D-l  lists  the  measured 
power  levels  and  the  percentage  power  level  referred  to  the  input  power. 

TABLE  D-l.  TRANSMISSION  AND  ABSORPTION  OF  A 
0.95-cm  LONG  PROUSTITE*  AT  1.06  pm 

Incident  Transmitted  Absorbed  Calculated  Re- 
Average  Power  Power  Power  flection  Loss 


mW 

percent 

mW  percent 

mW  percent 

percent 

Ordinary 

Wave 

153.5 

100 

77.4 

50.  5 

11.1 

7.24 

40.2 

Extraordinary 

160 

100 

82.6 

51.6 

9.6 

6.0 

38.2 

Wave 

*  Proustite  is  prepared  for  Type  II  upconversion  of  10. 6  pm  with  the  1. 06-pm 
pump. 
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As  evident  from  Table  D-l,  2  percent  of  incident  power  for  the 
ordinary  wave  case  and  4  percent  for  the  extraordinary  wave  case  are  un¬ 
accounted  for.  While  it  is  possible  that  this  missing  power  is  due  to  the 
measurement  error,  it  was  felt  that  the  scattering  inside  the  crystal  is 
responsible  for  the  bulk  of  these  losses.  The  accuracy  of  the  power  absorp¬ 
tion  measurement  is,  according  to  4.  Muray  of  Holobeam,  within  *5  percent 
of  actual. 

The  absorption  coefficient  of  the  material  was  computed  by  the 
following  equation: 

Pab-P0(,-c'Otl  (21) 

where 

P^b  ■  absorbed  power 

P  ■  power  input  to  the  material,  excluding  the  front  surface 
0  reflection  loss 

a  ■  absorption  coefficient 

l  m  material  length  in  cm. 

The  absorption  coefficient  for  the  ordinary  wave  was  found  to 

be  0.104  cm'1  and  that  for  the  extraordinary  wave  is  0.082  cm  The 

measured  absorption  coefficients  at  1.06  nm  are  one  order  of  magnitude 

3 

less  than  that  reported  at  694  nm  . 

During  the  absorption  measurement,  a  slight  reversible  tem¬ 
perature  dependent  absorption  was  noted.  Quantitative  measurements  of 
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the  temperature  dependence  were  not  performed  because  of  limitations  in 
the  instrument. 


Light  absorption  of  proustite  at  1. 06  pm  is  slightly  less  than 
-1  27 

the  0.45  cm  for  pyrargyrite  and  is  considerably  better  than  that  of 
ZnGePg^,  selinium,  and  CdGeASg^. 

The  temperature  dependence  of  the  light  absorption  and  the 
thermal  properties  of  the  proustite  should  be  measured  to  provide  a  better 
understanding  of  damage  mechanism.  It  is  also  recommended  that  a  more 
extensive  and  more  thorough  investigation  of  the  laser  damage  of  proustite 
be  done  under  controlled  conditions. 
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APPENDIX  E 

ANALYSIS  OF  HEAT  CONDUCTION  IN  LIGHT  ABSORBING  SOLIDS 


The  problem  of  heat  conduction  in  a  light  absorbing  uniform  iso¬ 
tropic  solid  is  considered.  A  nonlinear  material  such  as  proustite  used  for 
nonlinear  optics  application  is  generally  in  a  single  crystal  form.  Thermal 
transportation  in  the  single  crystal  is  not  isotropic,  but  the  tensor  properties 
of  thermal  conductivity  will  be  neglected  in  this  analysis  to  simplify  the 
problem.  It  is  also  assumed  that  light  absorption  and  propagation  are  iso¬ 
tropic  in  the  analysis. 

Figure  E-l  shows  schematically  the  problem  considered.  The 
material  is  cylindrical  with  its  side  firmly  in  contact  with  a  heat  sink.  The 
laser  beam  travels  along  the  cylindrical  axis. 


V//V////  ’ 


LASER  BEAM 


U„  +6(o) 


'  /U0  ♦ 6  (1 ) 


'///  U  //// ' 
0 


1(r.o)  e  ~az  (LASER  INTENSITY  INSIDE  THE  SOLID) 
U0+  6  -  AMBIENT  TEMPERATURE  AT  THE  END  FACE 
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FIGURE  E-l.  HEAT  CONDUCTION  GEOMETRY 


The  differential  equation  describes  heat  conduction  as: 


u  (r,  z,  t)  +  s  (r,  z,  t)  =  cp 


&u  (r,  z,  t) 


where: 


k  =  thermal  conductivity 


u  (r,  z,  t)  =  temperature  at  the  point  (r,  z)  and  at  time  t 


s  (r,  z,  t)  =  heat  source  within  the  medium 


where: 


c  =  specific  heat 


p  =  density 


r  =  transverse  coordinate 


z  =  axial  coordinate 


The  boundary  conditions  are  expressed  by 


u  (a,  z,  t)  =  U 


=  |-[u  (r,  o,  t)  -  9  (r,  o,  t)  -  UQ] 

z=o 


§1  =5  [“  fcM)  -  0  (r,t,t)  -  uo] 

Z  =  ^ 


l  =  axial  length  of  the  solid 


a  =  radius  of  the  cylindrical  solid 


•  ^  ^ - - 


(E-l) 


(E-2) 


(E-3) 


(E-3) 


•••  ■  •  ish1”  T 


Uq  =  temperature  of  heat  sink 

j 

h  =  coefficient  of  heat  exchange 
0  +Uq  =  temperature  of  medium  just  outside  of  the  end  faces 

The  boundary  condition  expressed  by  equation  E-2  indicates  heat  sinking  at 
the  side  wall.  Equations  E-3  and  E-4  indicate  heat  exchange  at  the  end  faces 
with  the  surrounding  medium. 

The  source  term  in  equation  E-l  is  given,  for  an  axially  propa¬ 
gating  light, 


s  (r,  z,  t)  =  iQ  (r,o,  t)  a  e  az  (E-5) 

where: 

i  (r,  o,  t)  =  light  intensity  at  the  input  end-face 
a  =  light  absorption  coefficient 

A  solution  to  equation  E-l  subjected  to  the  boundary  condition 
equations  E-2,  E-3,  and  E-4  can  be  separated  into  two  terms: 

u  (r,  z,  t)  =  v  (r,  z,t)  +  UQ  (E~6) 
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0 

0 

The  equation  and  the  boundary  conditions  for  v  then  become: 


fi 

.  2  by 

ktf  v+s=cpgj- 

(E-7) 

0 

v  (a,  z,  t)  =o 

(E-8) 

0 

by 

bz 

=  |  [v  (r,  o,  t)  -  0  (r,  o,  t)] 

(E-9) 

n 

z=o 

ii 

0 

_  5v 
bz 

=  |  [v  (r ,  l,  t)  -  0  (r,t,  t)] 
z= 

(E-10) 

11 

If  one  expresses  v  and  s  in  temporal  Fourier  integrals: 

n 

v  (r,  z,  t)  =  fv  (r,  z,  f)  e~l2ntt  df 

(E-ll) 

0 

s  (r,  Zj  t)  =  j S  (r,  z,  f)  e”i27rft  df 

(E-12) 

ii 

The  differential  equation  for  Fourier  component  V  becomes: 

ii 

ky2  V  +  S  =  -i2*fcp  V 

(E-13) 

i. 

The  boundary  condition  for  V  becomes: 

h 

V  (a,  z)  =  0 

(E-14) 

Eigenvalues  Xn  are  determined  by  the  boundary  cond.tion 


R  (a)  =0 
n 


The  solution  Rn  is  then 

R~  ^  =  Jo  ^yn  (E-22) 


where  y  is  the  nth  root  of  the  zeroth  order  Bessel  function  J  (x).  The 
n  ° 

source  function  S  can  be  expanded  in  terms  of  eigenfunctions  Rfl 


S  (r,  z,  f)  =  cl  e~az  ^ 
n=l 


H  (f)  R  (r) 
n  n 


(E-23) 


where 


IQ  (r,  o,  f)  Rn  (r)  rdr 


/3. 

R2  (r)  rdr 


The  differential  equation  for  ZR  (z)  is  obtained  by  substituting  equations  E-20 
and  E-  23  into  equation  E-13, 


d2  Z 


dz 


-  +  ( i2fffcpA  -  X2  ]  zn  =  -Hn  ae 


-az 


(E-24) 


The  boundary  conditions  for  (z)  are 


where 


d  Z 


n 


dz 


d  Z 


=  £-[Zn(°)-Tn(o)] 


z=o 


n 


dz 


■f  tZn  W-V*>] 


Z=^ 


T  ir,z)l  .  =[T  (z) |  ,  R  (r) 

v  *  '  z=o,l  n  '  /lz=o,-t  n  '  ' 


f 


T  (z)  |  ,  =- 

n  z=o,-t 


T  (r,  z)  Rn  (r)  rdr 


/  R*(r)rdr 
0 


z=o,l 


(E-25) 


(E-26) 


(E-27) 


The  solution  Z^  (z)  of  differential  equation  E-24  can  be  written  as 


Z  (z)  =  ^  (z)  •  e 

n  '  '  n  '  ' 


-az 


(E-28) 


Substituting  equation  E-28  into  equation  E-24,  we  obtain  the  differential 
equation  for  ^  (z) 


d  *  <W  „  on 
dz 


-H  a 
n 


(E-29) 


E-7 


where 


p  =  2fffcpA 


The  boundary  conditions  for  *  (z)  now  becomes 


d  * 


n 


d  z 


z=o 


i!" 

dz 


■I— 


The  solution  of  equation  E-29  is 


-II  a 

Vz)=_2 - 2"  2  ~*V 

or  -  y  /a  +  Ip 


"»»\b  > 

n 


where 


qnl  'f  4<4  vn/a2  -  3"2  •  i4p) 

0„2=f-l(4’'n/!,2-3“2-Wp) 


1/2 


1/2 


(E-30) 


(E-31) 


(E-32) 


(E-33) 


(E-34) 


Linear  simultaneous  equations  for  An  and  Bn,  which  obtained  by  substituting 
equation  E-32  into  equations  E-30  and  E-31,  are 


-H  o\a 
n  ,  k 

a  -  y  /a4  +  ip 
n 


ST„*» 


(E-35) 


+  B 

n 


4 


-H  or[  or  -  e^Jj- 
n  1  k 

<?  -  V^/a2  7  ip 


h  rrt 


+  KC  V4> 


(E-3G) 


It  must  be  remembered  that  the  function  *n  (z)  is  a  function  of  the  temporal 

frequency  f,  since  the  quantities  p,  H  ,  T  (o)  and  H  (<•)  are  all  functions  of 

n  n  n 

frequency. 

Solving  equation  E-35  and  E-3G  for  A  and  B  .  we  obtain 

n  n 


•  * 


-V* 


fc  ml 


A,  •  — 


h\ 

'  k  h 


H  • 


»  M  *  "  ' 

V"i 


•"Srt 


•t, 


V" 


CE-37) 


E-9 


where 


^-(^1-“-^!  (  (E-39) 


This  completes  formal  solution  of  heat  condition  in  light  absorbing 
solid  cylinder. 


r 
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